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£jj Abstract 

Consider a Hermitian matrix model under an external potential with spiked external source. 
When the external source is of rank one, we compute the limiting distribution of the largest 
eigenvalue for general, regular, analytic potential for all values of the external source. There is 
i— | a transitional phenomenon, which is universal for convex potentials. However, for non-convex 

potentials, new types of transition may occur. The higher rank external source is analyzed in 
the subsequent paper. 

g 1 Introduction and results 

^ 1.1 Introduction 

^-j- Fix an n x n Hermitian matrix A n and consider the following density function on the set H n of 

n x n Hermitian matrices: 

^ Pn(M) = ±r e -^(V(M)-A„M) (1) 

where Z n is the normalization constant. Here the 'external potential' V(x) is a real-valued function 
which decays fast enough as \x\ — > oo so that Z n is convergent. The matrix A n is called the external 
source: see e.g. [H], p2], [26], [27], [9], [TO], [2]. Note that the distribution of eigenvalues of M is 
unchanged if A n is replaced by UAnU^ 1 for any unitary matrix U . Since we are only concerned 
on eigenvalues of M , we assume without loss of generality that A n is a diagonal matrix. 



X 

A special case is when for all n, the external source has a fixed number m, called the rank of A 



of fixed non-zero eigenvalues. In this case, the sequence of probability spaces (H n ,Pn) is called a 
Hermitian matrix model with spiked external source, spiked source model for short. In this paper 
we only consider the case when m = 1. The higher rank case when m > 1 will be analyzed in the 
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upcoming companion paper. Throughout this paper, we assume that n > 1 and 

A n = diag(a, 0, • • • ,0), (2) 

n-l 

where a is a real number, independent of n. 

There are two important special cases. When V(x) = x 2 /2, the spiked source model is called the 
GUE spiked model. The density p n (M) is that of M = H + A n where H is an n x n GUE (Gaussian 
unitary ensemble) matrix. When V{x) = ((l+c)x — clogx)x(p,oo)( x )i c = (m—n)/n > 0, the spiked 
source model is the complex Wishart spiked model. In this case, setting £ := (1 — (1+c) A n ) _1 , the 
density p n (M) is that of M = Y^^XXW 2 where X is an n x m complex rectangular matrix with 
i.i.d. standard complex Gaussian entriesj^] For these two cases, the limit of the largest eigenvalue 
Cmax(^) of M was studied in great detail in [3] and [23]. An important feature is the following phase 
transition phenomenon. Let e denote the right-end point of the limiting empirical distribution of 
the eigenvalues of the Hermitian matrix model with no external source (see ([10} below) It was 



shown in both the GUE and the complex Wishart spiked models that as n — > oo, with probability 
1, 

U«(n) -> < u (3) 

for some continuous, increasing function Xq{o) in a G (^V'(e),oo) satisfying lim^iy,^^ Xq{q) = e. 




Moreover, there exists (3 > (see ( 22 ) below) such that for each T £ 

if o 

if a = \V'(e) 
and there exists 7(a) such that for each T € M, 

Pn((£max(n) - x (a)) 7 (a)n 1 / 2 < T) — >■ G(T), if a > \v\e). (5) 

Here the function G(T) = /j 7 ^ e~2? is the cumulative distribution function of the standard 
normal distribution, and Fq and y/F\ are the GUE and GOE Tracy- Widom distribution functions, 



respectively. They are defined in (16) and (|20j) below, respectively. A limit theorem was also proven 
for the double scaling case when a = ^V'(e) + -573. 

The purpose of this paper is to extend the results ([3])-([5]) to the spiked source model with general 
potential V. It turns out that if V(x) is convex in the interval x G (e, 00), then all of — ([5]) still 
hold. Especially, the 'critical value' of o is again given by ^V'(e). However, if V is not convex in 
(e, 00), new features may occur. Two key new features are the followings. 



1 For the complex Wishart spiked model, V is not real analytic at x — 0. Throughout this paper, we only consider 
V which is real analytic in the whole line. However, the method can be generalized to the Wishart-type potentials 
in a straightforward way. 

2 The limiting empirical distribution in the spiked source model is the same as the Hermitian model with no 
external source. 



2 



1.2 



and Theorem 



1.2 



For 



The critical value of a may be smaller than hV'(e). See Lemma 
such a case, when a equals this critical value, ( ma); (ii) does not converge with probability 
1. Instead it converges to two or more values, each with non-zero probability. In this case, 
the fluctuation of £ max (n) is generically Fq at the smallest limiting value and G at the larger 
limiting values. See Theorem |1.3| 



• There may be a discrete set of 'secondary critical values' of a, which are greater than the 
critical value. If a is at a secondary critical value, then £ max (^) converge to two or more 
values, each with non-zero probability. In this case, the fluctuation of ^ mH (n) is generically 



G at each of the limiting values. See Theorem 1.4 



The exact assumptions on the potential V is given in Subsection 1.2. The universality result for 



convex potentials is in Subsection |1.3| In Subsection 1.4 we define the critical and the secondary 
critical values for non-convex potentials. The limit laws for the non-convex potentials are given in 
Subsection 11.51 



While we were preparing for this paper and the companion paper for the higher rank case, 
we learned that M. Bertola, R. Buckingham, S. Y. Lee and V. Pierce were also working on the 
spiked source models (see [8] for the first part of their work). While we focus, especially in the 
second paper, on the limit laws when ai, • • • , a m are distinct, Bertola, Buckingham, Lee and Pierce 
focus on the case when a\ = ■ ■ ■ = a m and m — > oo slower than n. Also we use the asymptotics of 
usual orthogonal polynomials but Bertola, Buckingham, Lee and Pierce use asymptotics of multiple 
orthogonal polynomials via Riemann-Hilbert problem of size larger than 2. 

Before closing this subsection, we mention that the spiked real symmetric matrix model is much 
more difficult. Even for the GOE and the real Wishart case, the limiting distribution at the critical 
value is not yet known. For the quaternionic case, the limiting distribution is obtained when the 
rank m = 1 (see |25j for the Wishart model; Gaussian model is also similar). 

We also mention that there are several results for the spiked Wigner ensembles and spiked 
sample covariance matrices. See, for example, [S], [22], [T5], [13], [2T], [7J and [BJ. 

1.2 Assumptions on external potential V. 

Throughout this paper, we assume the following three conditions on V: 

V{x) is real analytic in K, (6) 
V(x) 



, , ( ■> +oo as \x\ — > oo, (7) 
\J x 1 + 1 

V is 'regular'. (8) 

The second condition is to ensure the convergence of the density function: compare this with the 
condition on V in |17j . The third condition on being 'regular' is a technical condition as defined in 
[17] . We need a few definitions to state it. 
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First, recall the equilibrium measure and the so-called g-function. General references are 
and [15J. For a given potential V, the empirical distribution of the eigenvalues of the matrix 
model with no external source converges to the associated equilibrium measure /i. The equilibrium 
measure is characterized by a certain variational problem. If V is real analytic, [i is supported on 
a finite union of intervals, 

JV 

J = U^^-kOj with b < a\ < b\ < ■ ■ ■ < a N+1 , (9) 
3=0 

for some N > 0. We denote the right-most edge of the support by 

e := ojv+l- (10) 

On J, d[i has the form d[i = ty(x)dx, 

1 N 

«r(a?) = R 1 / 2 (x)h(x), where R(z) = JJ(z - bj)(z - a j+1 ). (11) 

The function R(z) 1 ^ 2 is defined to be analytic in C \ J and satisfy R(z) l l 2 ~ z N as z — > oo. The 
notation R^ 2 (x) for x G J denotes the limit of R 1 ^ 2 (z), z 6 C+, as z — > x from above. The function 
h(x) is real analytic in M and is given by [17} Formula (3.18)]. 

The equilibrium measure, d/j,(x) = ^f(x)dx, is characterized by the following conditions: there 
is a constant (called the Robin constant) £ such that 

2 J log\x - s\^(s)ds - V{x) =£ for x G J, (12) 

2 J log\x - s\^(s)ds - V(x) <t for x G M \ J. (13) 
The so-called g-function is defined by 

g(z) := j log (2 - s)ty(s)ds, for 2 G C \ (-00, e). (14) 

The potential V is said to be regular (see [17]) if 
• h(x) 7^ for x G J, 



• the inequality in (13) is strict. 

The first condition implies that the function ^(x) > for all x G J, and also that ^(x) vanishes 
like a square-root at each end of the interval of the support. This in turn implies, in particular, 
that for the model with A n = 0, the largest eigenvalue has the limiting distribution given by Fq 
(see e.g. [J3] for the non- varying weight; varying weight case is similar using the analysis of |17j.) 
Note that the second condition restricted to the domain x > e implies that 

2g(x) - V(x) <£, x > e. (15) 

We will use this fact later. 
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Figure 1: The contour from ooe 57 ™/ 6 to ooe 7 ™/ 6 . 

1.3 Statement of results: convex potentials 

Let Fq(T) be the GUE Tracy- Widom distribution defined by 

F (T) := det(l - X[T,oo)-^Air y X[T,oo)), (16) 

where X[T,oo) denotes the projection operator on [T, oo), and Kj^ y is the Airy operator defined by 
the kernel 

K-Ahy{X,y)= • (17) 

x — y 

Here Ai is the Airy function. 

For qGR, define the function 

C a (0 e*****-^-, (18) 

ZTT J a + iz 

where the contour is from ooe 57 ™/ 6 to ooe 7 ™/ 6 and the pole z = —ia lies above the contour in the 
complex plane: see Figure [TJ 
Define 

F 1 (T;a) := F {T) ■ (l - ((1 - XlT^^AiryXlT.oo))" 1 ^, Ai) [TiCX)) Y (19) 

where (/, g) e denotes the real inner product over E, j E f(x)g(x)dx. (See [31 Definition 1.3].) When 
a = 0, 

Fi(T) := Fx(T; 0) (20) 

equals the square of the GOE Tracy- Widom distribution (see [31 Formula (24)]). 

Fix a potential V satisfying the assumptions ([6])-([8]). In the companion paper on the higher 
rank case, we need to consider the spiked source model of rank 1 whose density function is same 
as in ([!]) but with the change that the matrix M is now of size n — j + 1 and A ra is replaced by 
A n _ J+ i for fixed j: 

1 _ e -nTr(V(Af)-A n _ J+1 M) ) M € U n -j+l. (21) 



Z n -j+l,n 

Note that the factor n in front of the potential is unchanged. For a subset Bel, let P n _j+i jn (a; E) 
denote the gap probability that there are no eigenvalues of M in the set E, where a represents the 
unique non-zero eigenvalue of A n _ J+ i. Hence P n _j + i >n (a; [t, oo)) is the probability that the largest 
eigenvalue of M is less than t. 
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Let e be as in (10). Set (recall (11)) 

vr^(x) \ 2/3 



lim , 

arte ye — X 



He) 



2/3 



R(z) 



z — e 



1/3 



(22) 



so that ^(x) 



£3/2 



\/e — x for x f e. For T£l, define the intervals 



and 



I T :-- 



e + 



/3n 2 / 3 



, oo 



T 



, oo 



(23) 



(24) 



^(V"(x*)-%"{x*))n 

for x* > e, if V"(x*) — g"(x*) > 0. Note that if V(x) is convex in x > e, then V"(x) —g"(x) > for 
all x > e. For later reference we note that V"(x*) — g"(x#) = — G"(x*) in terms of the notation (32 ) 
that is defined below. 

The following is the first main result of this paper. Let V(x) be a potential that is convex in 
x G (e, oo). For a > |V'(e), let xo(a) be the unique maximizer of the function g(x) — V(x) + ax 
in x G (e, oo). Such a maximizer exists since g(x) — V(x) in x G (e, oo) is strictly concave and 
g'(e) — F'(e) + a = — ^V'(e) + a > (see (p0|)) and g'(x) — F'(x) + a < for all large enough x. 



This xo(a) is same as in Lemma |l.3[ 

Theorem 1.1 (convex potential). Let V(x) be a potential that is convex in x G (e, oo). Set 

a c := \V(e). 

The following holds for each T G K as n —> oo anc? j = O(l). 
(a,) For a < a C; 



(25) 



(7>j For 



lim P n -j+l,n (a; /J) = F (T). 

/3a 



a = a c + 



1/3' 



where a is in a compact subset of M., we have 

lim P n _j + i n (a; /J) = Fi(T; -a). 

(cj For a > a c , 

lim P n _ i+ i, n (a; jj(x (a))) = G(T). 



(26) 
(27) 

(28) 
(29) 



Hence the transition phenomenon is universal for convex potentials. The next two subsections 
are about non-convex potentials 
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1.4 Critical value and secondary critical values 

In this Subsection, we define critical values and the secondary critical values of a. 



By definition (14), g(x) is real analytic in (e,oo), is continuously differentiable in [e, oo) and 
satisfies 

g'(x) > 0, g"(x) < 0, forxG(e,oo), 

V(e) + £ V'{e) , (30) 

g(e) = - , g( e ) = — — , limg(x) = 0. 

i 2 x— >oo 

Definition 1.1. For a G (0, ^V'(e)), define c = c(a) as the unique point in (e, oo) satisfying 

g'(c(a)) = a. (31) 

For a > 7jV'(e), define c(a) := e. 

Note that c(a) decreases strictly in a £ (0, V (e)/2) and continuous in a G (0,oo). 
Define two auxiliary functions 

G(z) = G(z;a):=g(z)-F(z)+az, 
H(z) = H(z; a) := -g(z) + az + t 

for z£C \ (— oo,e). Observe the following Lemma. The proof follows straightforwardly from the 
definition of g, the variational condition (12), the assumption (j7|) on V and (15). We omit the 



details. 

Lemma 1.1. Let a > 0. We have the following properties. 

(a) H(x) is a convex function in x G [e,oo) with the unique minimum attained at x = c(a). 

(b) H(x) > G(x) for all x G (e, oo). 

(c) H(e) = G(e) = -\V(e) + ae + \l. 

(d) limH'(x) = limG'(x) = a - -V'(e). 

x\.e x\.e 2 

(e) As x — > +oo ; H(x) — > +oo, H(x)/x — )• a, G(x) — > — oo anc? G(x)/x — > — oo. 
See Figures [2] [3] and [4] for a few examples of the graphs of G and H. 

Define the set 

Ay '■= {a G (0,oo) I there exists x G (c(a), oo)such that G(x;a) > H(c(a);a)}. (33) 
Definition 1.2. The critical value for the spiked source model with potential V is defined as 

a c :=inf^ly. (34) 
Lemma 1.2. We have the following properties. 
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c(o) = \V'{e) x (a) 
a > a c 



Figure 2: Schematic graphs of the functions H(x;a) and G(x;a) for a potential V such that 
a c = ^V'(e), assuming that a $ Jy. 

(a) (±V'(e),oo) C -4y #ence a c < ^'(e). 

^ The set Ay is an open, semi-infinite interval. Hence Ay = (a c ,oo). 

(c) a c > 0. 

(d) For < a < a c , we aaue G(x;a) < H(c(a);a) /or aZZ cc G (c(a), oo). 

(e,) // the potential V(x) is convex for x > e, then a c = |V'(e), and G(x; a c ) < H(e; a c ) /or all 
x > e. (Note that c(a c ) = e and G(e; a c ) = H(e, ac) .) 

(f) If the potential V is such that a c < | V'(e), £/ien G(x; a c ) < H(c(a c ); a c ) /or aZZ x G (c(a), oo) ; 
and f/ie equality is attained at least at one point. 

Proof. |(a)| Let a G (^V'(e),oo). Since lim^e G'(x) = a — ^V'(e) > 0, there is x > e such that 
G(x) > G(e) = H(e). Thus a G >4y. 
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Figure 3: Schematic graphs of the functions H(x;a) and G(x;a) for a potential V such that 
a c < ^V'(e), assuming that a g" Jy. 



|(b)| The continuity of G and H in a implies that Ay is an open set. Now we show that A is a 
semi- infinite interval. Suppose that a G Ay and a < hV(e). Let x G (c(a), oo) be the point 
such that G(x;a) > H(c(a);a). Let a' G (a, ^V'(e)]. From Definition 1.1 of c(o), we see that 
c(a') < c(o), and hence x G (c(a'), oo). Moreover, 

G(x;a')-U(c(a');a') = [G(x-,a)-U(c(a)-,a)} + [li{c(a);a')-U(c(a')-,a')] + [{a' -a)(x-c{a))} 

(35) 

is strictly positive since each term in bracket is strictly positive. Thus a' G Ay, and this, 
together with (a), implies that Ay is a semi- infinite interval. 



We have G(x) - H(c(a)) = -V(x) + a(x - c(o)) + g(x) + g(c(a)) -t < -V(x) + ax + 2g(x) -I. 
This tends to — oo as x — > +oo due to the growth condition ([7]) on V. Also c(a) —> +oo as 
a I 0. Therefore when a is close to 0, G(x) — H(c(a)) < for x > c(a). Hence, a .Ay if a 
is small enough. 

|(d)| Let < a < a c . Suppose that there is x G (c(a), oo) such that G(x;a) = H(c(a);a). For 
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Figure 4: Schematic graphs of functions H(x; a) and G(x; a) when a = a c G .Ty an d a c > a G 



any a' G (a,a c ), we have c(a) > c(a') since a < a' < ^'(e). Thus we find from (35) that 
G(x; a') — H(c(a'); a') > 0. This implies that a' G Ay which is a contradiction. 



(e) Let < a < ^V'(e). We will show that a ^ A. Since V is convex, G(x) is concave in 
x G (e, oo). As G'(e;a) < 0, this implies that G(x) is decreasing in x G (e, oo). Thus for 
x G (c(a), oo), G(x) < G(c(a)) < H(c(a)). Hence a ^ „4y. When a = ^V(e), a similar 
argument implies that G(x) < H(x) for all x > e. 

|(f)| This follows from the continuity of G and H in a and the fact that a c = inf Ay ■ 

□ 

See typical graphs of G and H for a c = ^V(e) in Figure [2J and typical graphs of G and H for 
< a c < \V'{e) in Figure g 

Remark 1.1. When V is non-convex, there may exist x > e such that G(x;a c ) = H(e,a c ) even if 



a, 



§V(e). 



By Definition |1.2| of a c , when a > a c , G(x;a) > H(c(a);a) for some x > c(a). The point 
x at which G(x; a) attains its maximum plays an important role. Indeed, we will show in the 
below that if the maximum is attained at a unique point, then ( max (n) converges to this point (see 



Theorem 1.2). However, it may happen that for some a's, the function G(x; a) attain its maximum 



at more than one point. Let 



c (a) := max G(x;a), (36) 

ic£[c(a), 00) 



and define 

J v ~{a G [a c ,oo)| G max (fl) is attained at more than one point}. (37) 

This set is discrete since G(x; a) is analytic in both x and a. Note that when V is convex, a c g' Jy 
from Lemma 1.2|[e)| For a non-convex V, as indicated in Remark a c may or may not be in Jy 



no matter if a c = ^V(e). See typical graphs of G and H for a G Jy in Figure |4j 
We have the following Lemma. 
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Lemma 1.3. (a) For a G [a c ,oo) such that a G" Jy, let xo(a) be the unique point in [c(o),oo) at 
which G(x;a) attains its maximum. Then xo(a) is a continuous, strictly increasing function 
in a G [a c , oo) \ Jy. 



(b) If ao G Jy and ao > a c , then 



limxo(a) < limxo(a). (38) 

a'fao aXdQ 



Note that if a G Jy satisfies a > a c or a = a c < |V (e), then there exist points x\(a) < X2(a) < 
■ ■ ■ < x r (a) in (c(a), oo), for some r > 2, such that 

Gmax(a) = G(xi(a); a) = • • • = G(x r (a); a). (39) 

On the other hand, if V is a potential such that a c = ^V'(e) and a c G J7y, then there exist, for 
some r > 1, xi(a c ) < X2(sl c ) < • • • < x r (a. c ) in (e, oo) such that 

G(e; a) = G max (a) = G(xi(a c ); a c ) = • • • = G(x r (a c ); a c ). (40) 



Proof of Lemma 1.3. The continuity of xo(a) for a Jy is a direct consequence of the continuity 
of G(x; a) in both x and a. Let a c < ai < and ai, 02 i7y. If we assume xo(ai) > xq{o2), then 
since G(rro(ai); ai) > G(xo(a2); ai), we have 

G(x (ai); 02) = G(x (ai); «i) + (a 2 - ai)x (ai) > G(x (a 2 ); a\) + (a 2 - ai)x (a 2 ) = G(x (a 2 ); a 2 ). 

(41) 

This is contradictory to the assumption that a?o( a 2) is the maximizer of G(x;a2). Thus xo(a\) < 
x (a 2 ). 

If ao G Jy and ao > a c , then G(x;ao) attains its maximum in [c(ao),oo) at a;i(ao), . . . ,x r (ao) 
for some r > 2 as in (39 ). It is easy to check from the continuity of G(x; a) in a that lim a | ao xo(a) = 



xi(a ) and lim a ^ ao x {a) = x r (a ). □ 

Definition 1.3. The secondary critical values for the spiked model are defined as the points a G Jy 
such that a > a c . 

Remark 1.2. For a potential V such that V(x) is convex for x > e, Jy = since G'(x;a) is a 
decreasing function in x > e. Hence there is no secondary critical value. 

1.5 Statement of results: non-convex potentials 

Let V{x) be a potential satisfying the conditions ([6])-([8]). Let /J and Jj(x*) be the intervals 



defined in (23) and \24m, respectively. 



Theorem 1.2 (away from critical values). The following holds for each T G M. as n — > 00 anc? 
3 = 0(1). 

(a) For a < a c , 

lim Pn-i+l,n (a; /J) = Fo(T). (42) 
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(b) For a > a c such that a £ Jy, if G" (xq(o)) ^ 0, then 

lim P n -;+l,n (a; (a))) = G(T). (43) 

When a is at or near the critical value a c , we have the following result. The case when a = a c 
is attained by setting a = 0. 

Theorem 1.3 (at or near the critical value). We have the following for each T6l. 
(a) Suppose that V is a potential such that a c = hV'(e) and a c ^ Jy. Then for 

a = a c + — , (44) 
where a is in a compact subset ofM, we have 

lim P n _ i+ i, n (a; /J) = F\(T\ -a). (45) 

f&j iet V be a potential such that a c < ^V'(e). //a c ^ J7i/ anc ^ G"(xo(a c ); a c ) / 0, then for 

a = a c + -, (46) 
n 

where a is in a compact subset ofM., we have, as n — ^ oo ; 

P„_ i+ i, n (a; /J) = p j>n (a)F (T) + o(l) (47) 

and 

P„- i+lj „ (a; jj(x (a c ))) = Pi , n (a) + (1 - Pj> )(a)G(T) + o(l), (48) 

where the constant pj >n (a) S (0, 1) is defined by ( |145 ). ^4s a function of a, pj, n {a) is decreasing 
and satisfies pj )U (a) — > as a —> oo and pj jn (a) — > 1 as a — > — oo /or eac/i ,/ixed n. Also for 
each fixed a, pj n (a) lies in a compact subset of (0, 1) for all large n. 

Remark 1.3. When the potential V(x) is convex for x > e, then a c = |V(e), J7y = (see Re- 
mark 1.2 ) and G"(x) < for all x > e. Hence Theorem 1 1 . 2 1 and Theorem 1.3fta) imply Theorem |l.l| 

For a at or near the secondary critical values Jy \ {a c }, we have the following result. 

Theorem 1.4 (at or near the secondary critical values). Let V be a potential such that Jy ^ 0. 
Let ao G Jy \ {a c } be a secondary critical point. LfG(x; clq) attains its maximum G max (ao) at two 
points xi (ao) < X2(ao) in (c(a), oo) and if G"(xi(ao); ao) ^ and G"(x2(ao); ao) ^ 0, then for 

a = a + -, (49) 
n 

where a is in a compact subset ofM, we have 

Pn- J+ i,n (a; Jj(*i(a ))) = pS(«)G?(T) + o(l) (50) 
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and 



j+1 , n (a; Jj(x 2 (a ))) = p^(a) +^(a)G(T) + o(l) 



where p*p n (a) and G (0,1) are defined in (164) and (158), and p~^(a) +pV^(a 

function of a, pf\(a) is decreasing and satisfies p^Aa) — ^ as a — )■ oo andp^(a) — ^ 1 as a 



.(2) 



(1)/ 



(51) 



1. vis a 



-oo 



/or eac/i /lied n. Also for each fixed a, p^{a) is in a compact subset of (0, 1) independent of n. 

The above three theorems describe the 'generic' cases. The next part describes the three 
'exceptional cases'. 

As the first exceptional case, suppose that in Theorem 1.4 the maximum of G is attained at 
more than two points. Let x\(a) < x 2 (a) < • • • < x r (a) be these maximizers. If G"(xj(ao); a>o) ^ 
for all k = 1, • • • , r, then we have, for each k = 1, • • • , r, 



fc-i 



n-j+l, 



,„ (a; J n T (x fc (a ))) = £>$(a) +p^( a )G(T) + o(l) 



(*)/ 



(52) 



i=l 



for some p^( a ) G (0, 1) such that j>$(a) + • • • + P^(a) = 1. Explicitly, ^>(«) ■= Al ( a )+-+A r {a) 
where Ai(a) is defined in (158). The situation when a c 6 J7y in Theorem 1.3;[b) is similar. In this 
case, the maximum of G(x; a c ) in (c(a c ), oo) is attained at ^i(a c ) < X2(a c ) < • • • < x r (a c ) for some 
r > 2 (see (39)). Assume that G"(xj(a c ); a c ) / for all i = 1, • • • , r. Then with Cj(a), i = 1, • • • , r, 

1, • • • , r, where 

>— -P$(°). 



,(*) 



defined by (144) withxo(a c ) replaced by 2j(a c ), setpf^ (a) : 



Cj(a) 



C +Ci(a)+-+C r (a)> 

Co is defined by (143). Then (47) holds with pj >n (a) replaced by p^(a) := 1— p^}\(ot 
The limit in (48) is replaced by 

fc-l 

P n _ i+1 , n (a; J n r (x fe (a ))) = J>$(a) + P$(a)G(T) + o(l) 



i=0 



(53) 



for fc = 1, ■ ■ ■ , r. 

The second exceptional case is when a c G J7y in Theorem 1.3, case 1. This case is given in the 
following Theorem. In this case, there are two natural scalings in a. 

Theorem 1.5. Let V be a potential such that a c = ^V'(e). Suppose that a c G Jy . Assume that 
G(x;a c ) attains its maximum at the unique point xo(a c ) G (c(a c ),oo) and G"(xo(a c ); a c ) ^ 0. 
Then the following holds. 



(a) For 



a c + 



/3a 
^73' 



n 



where a is in a compact subset of (— oo,0), we have 

lim F n (a-ll) =F 1 (T;-a). 



(54) 



(55) 
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(b) For 



a c + 



a 



n 



where a' is in a compact subset ofM, we have 

Pn-;+l,n («; In) = Pj,n( a ') F l( T ) + °(!) 

where *i(T) = Fi(T;0), and 

V n - j+ i, n (a; Jj(x (a c ))) = p j)n (a') + (1 - p j>n (ct))G{T) 



(56) 
(57) 
(58) 



where pj iTl (a') is defined in ( |280 ). ^4s a function of a, pj tn (a!) is decreasing and satisfies 
Pj,n(oi') —7-0 as a' —7- oo and pj tn (a') — > 1 as a' —> — oo /or eac/i /ixed n. ^4teo /or eac/i /ixed 
a, Pj jn (a) is in a compact subset of (0,1) independent of n. 



If the maximum of G(x;a c ) is attained at more than one point, then (58) should be changed 
in a natural way as in ( |53| ). 

The third exceptional case is when the double derivative of G(x; a) vanishes at its maximizers. 
Then the function G(x) is replaced by its higher analogue, and the scalings in the interval and a 
are also changed accordingly. Concretely, in Theorem l.^[b)[ if G"(xo{a)) = 0, then since xo(a) 
is the maximum point, there exists k > 1 such that G^ 2fe -'(xo(a)) < and G^(a;o(a)) = for all 
j = 1,2, ■ ■ ■ ,2k — 1. Then (43) is changed to 



lim ] 

n— >oo 



j+i,n («; J n (x Q (a);k) 



I- 



e ax 



'dx 



(59) 



where the interval </J(x*; k) is defined by 



X* ~\~ 



' n{V^ 2k \x,)) - ^ 2k \x,)) 



(2*)! 



T, oo 



(60) 



In Theorem 



l.c 



b) 



Theorem 



0, (i = in Theorems 1.3j[b) 



1.4 



and Theorem 



1.5 



a is scaled as a = «o + f • When G"(xi(ao); ao) 



and 1.5, and i = 1,2 in Theorem 1.4,) then this scaling also needs to 



be changed. For example, Theorem |1.4| is changed to the following theorem. 

Theorem 1.6. Let ao E Jy \ {a c } be a secondary critical value. Assume that G(x;ao) attains its 
maximum G max (ao) at two points xi(ao) < £2(00) in (c(a), 00). Suppose that G"(xi(ao); ao) 7^ 0, 
and for some k > 1, and suppose that G ( - 2k \x2(ao); ao) 7^ and GW(x2(ao); «o) = for all 
i = 1, . . . ,2k — 1. Then for 



a = ao 



logra a 

/ + -; 

n n 



where q 



2k 



x 2 (a ) - xi(a ) 
where a is in a compact subset o/M, we have 

P n _ j+1 , n (a; Jj0n(ao))) = p£(a)G(T) + o(l), 



(61) 



(62) 
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and 



n-j+i,n (a; J 1 n (x 2 (a ); k) J = pfl(a) + p)%(a) 



S 2 ), 



J-oo e 



x2k dx 



T- 



'dx 



(63) 



where the interval J* (x; k) is defined by (60), p^(a) and pV^(ai) are defined in (182), andpV^(a)-|- 



( 2 )/ 



Pjn( a ) = 1- ^ S a function of a, p^(a) is decreasing and satisfies p^(a) 
Pjn( a ) ~ * as a — ^ — 00 ^ or eac ^ /^ed n. /or eac/i /ixed a, pj-^a) is in a compact subset of 



as a -> oo and 



(0,1) independent of n. 



The changes needed for Theorem 1.3|[b) and Theorem 1.5 are analogous. Also it may happen 
that the two or more of the exceptional cases occur simultaneously. Then one needs simply combine 
the results together in a straightforward way, and we skip the details. 



We remark that if the support J of the equil ibrium measure is o f o ne in terval i.e. N 

and 

, * * J '" ~ : 11 1 II— 

follows from Remark 



6.1 



1.6 



(see ©), 
do not depend on n. This 



re definition of these probabilities. When 



on M.j, n , -Mj,n and Bj >n (e) and t : 
N > 0, the dependence of these probabilities on n is from the theta function formula of ^Aj^ n , ]\A^j 
and Bj >n (e) in Section [61 and is in a quasi-periodic way. 



We also remark that one can obtain the convergence in probability 1 as in ^ from the above 
theorems together with the fact that all of the limiting distributions decay rapidly at the tails. 

An explicit example of a potential such that a c < ^V'(e) can be constructed as follows. We use 
the potential defined in [18} Formula (4.14)] (we change the original notation e into e here): 



V, 



1 



[x 



l + ee V 4 



1 



e + e 



-x° + 



ee 



-x l + 2(e + e)x 



(64) 



where e is a very small positive number, e is any number > 2 and e is determined by e from the 
condition that 



J (x — e) (x — e) y/a 



2 - Adx = 0. 



(65) 



From results in [18} Section 4], it is known that Vg i£ is a regular potential with the support of the 
equilibrium measure given by [—2 + 0(e), 2 + 0(e)\. For all x > e = 2 + 0(e), ( 15 ) holds. However, 
at x = e, 

2g(e) - 7(e) - t = -E(e) (66) 

for some E(e) satisfying E(e) = 0(e) and E(e) > 0. Hence for any a, G(e; a) — H(e; a) = O(e). Now 
there exists a E (0, ^7'(e)) such that c(a) S (e, e) since the minimizer c(a) of H(x; a) is continuous 
in a G (0, oo), decreases strictly in a £ (0, |7'(e)), lim^o c(a) = +oo and c(hV(e)) = e (see the 
sentence after Definition [Li]) . Since H(c(a); a) < H(e; a), we have G(e; a) > H(c(a); a) if e is small 
enough. Then a S .Ay and a c < 57 g ' e (e) for each fixed e > 2 if e > is small enough. 

The paper is organized as follows. The outline of the proof of theorems is given in Section [2j 
The results on the orthonormal polynomials and the kernel K n _j^ n are summarized in Section |6j 
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The proofs of the theorems are given in Sections |3j |4| and |5j We consider three cases, a c < ^V'(e), 
a c > iy'(e) and a c = |V(e) separately. Throughout this paper we only consider a > 0. The 
a < case is discussed briefly in the end of Section [2} 
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2 Outline of the proof 

Let 

Pj(x; n) = 7j(n)j; : ' + • • • (67) 

be the orthonormal polynomial of degree j with respect to the weig ht e~ nV( - x \ Here take jj(n) > 
to make pj unique. Set 

i)j(x;n) := pj(x;n)e~% y ( x \ fj{x; n) := pj(x; n)e~ nV ^ x \ (68) 

Let 

Kj, n {x, y) := > ifri(x; n)ipi(x; n) = J - J - J - J - (69 

~^ lj\ n ) x-y 

be the Christoffel-Darboux kernel. Define the constant 

r>;n) := ( e «(«*-v(«)/2) > ^. {x . „ ))r> (70) 

and the function 

j>i(x;a;n) :=^—, ( e n{ax ~ V{x)/2) - £ ^(a; n)^(x; n)) 

r ^ a;n ) V U ) (71) 



1 ( n(ax-y(,)/2) _ /" # ( x>y ) e »(«y-v(i0) d 



r,(a; 

It is easy to check that ipj{x; a; n) is characterized by the orthonormality conditions 

(tjjj(x;a;n),ip k (x;n)) R = 5 jk , for k = Q,l,...,j 1 (72) 

in the vector space spanned by {e _ a v ( x \ xe~z v<yX \ • • • , xi~ 1 e~2 V ( x \ e n(ax-V(x)/2)y_ j n ^ e mu ltiple 
orthogonal polynomial theory, ipj (x; a; n)e nV ^' 2 is the multiple orthogonal polynomial of type I 
with potentials e~ nV<yX ^ and e ax , see e.g. |10j . From this follows the well-definedness of , ijjj(x;a;n), 
i.e., Tj(a; n) ^ 0. 

We sometimes drop the dependence on n or a in ipj(x;n), Tj(a;n) and ipj(x;a;n) and write 
iftj(x), Tj(a) and ^-(a;) for convenience. 



16 



The starting point of our analysis is the following basic result in the theory of Hermitian matrix 
model with external source, specialized to the spiked source model of rank l(see e.g. [26], |10j): for 

any Eel, 

P„-j+i,n(a; E) = det (l - XEK n -j+i,nXE) , (73) 

where 



j+l,n 



n—j,n 



+ lfj n -j ® 1p n - 



Here xe denotes the projection operator on E and Kj <n is the operator on L 2 (M 
the kernel (|69|). Note that the only term in (74) that depends on a is ijj n -j(x;a;n) 



(74) 

defined by 
The kernel 



K n -j,n{x,y) is precisely the reproducing kernel in the Hermitian random matrix model of size 
n — j with weight e~ nV with no external source. Hence for the rank 1 spiked Hermitian model, the 
reproducing kernel is a rank 1 perturbation of K n _j^ n . For the higher rank spiked Hermitian model, 
the reproducing kernel is a rank r perturbation of K n _j tn which will be studied in the subsequent 
paper. 



For the asymptotic result for this paper, we need asymptotics of K, 



n—j,ni 



tpn-i and V>n- 



The 



asymptotics of orthogonal polynomials and the Christoffel-Darboux kernel with respect to a varying 
weight have been studied extensively. Most notably precise strong asymptotics were obtained for a 
general class of potentials using the Deift-Zhou steepest-descent method for the associated Riemann- 
Hilbert problem (see [T7] and [H]). We use the results of p2] extensively. However, in [17], only 
the case of j = 0, 1 are stated explicitly. For more general j = O(l), the same analysis of p2] can 
be carried out after a few changes. These have been studied in various other papers (see e.g. [I] 
for the discrete orthogonal polynomials case). We summarize the necessary changes and state the 
explicit asymptotic formulas in Section [6] below. 

The main part of this paper is the asymptotic analysis of ip n -j(x;a;n) defined in (71). For 
this purpose, we first evaluate the asymptotics of T n ^j(a; n). This can be achieved in principle by 
plugging in the asymptotics of ip n -j(x;n) in the definition (70) of T n _j(a;n), and then evaluating 
the inner product asymptotically. However, the oscillatory nature of ip n ~j(x; n) in the support of the 
equilibrium measure makes it cumbersome to evaluate the inner product in this way. Instead we re- 
express T n -j(a; n) in terms of a sum of integrals involving both ^ n -j(x; n) and its Cauchy transform 
(see (85)). This removes the oscillation and the asymptotic analysis becomes more straightforward. 
A similar trick is also used in evaluating J R K n _j in (x, y)e n ^ ay ~ v ^dy (see Lemma 3.1 ). The analysis 
is divided into several cases depending on the location of the critical points. Each of these cases 
correspond to the theorems in Introduction. 

For the set E = or E = J^(x^) of interest in each theorem, we can show that ijj n -j(x; a;n) 
is in L 2 (E) for each n, and hence XE^n-j <8> ^n-jXE in (74) is a trace class operator. However, 
the L 2 -norm of ip n -j is not uniformly bounded in n nor there is a simple estimate on the L 2 norm. 
For example, we will find that in (124) that tp n -j(x; a;n) = 0(^/ne n( - G ^ +n ^~ 2H ^^ 2 ) for x > c 
when a < min{a c , ^'(e)}. The function G(x) + H(x) - 2H(c) = -V{x) + 2a(x - c) + 2g(c) - I 
tends to — oo as x — > +oo, but it may be positive for some value x > c. This implies that we 
do not have a good trace norm of XE^n-j <8> tp n -jXE = llV ; n-j||L2( £ ;)||'0 n _j|| L 2( £ ;- ) , and we cannot 
compare the size of XEK n -j,nXE and XEipn-j <S> tpn-jXE- However, the rapid decay of the operator 
(1 — XEK n -j,nXE)~ l can be used to control the estimates. We proceed as follows. 
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Prom (73) and (74), 



P n _ i+ i in (a; E) = det (l - XEKn-j+l,nXE) 

= det (1 - XE~K n -^ n XE) ■ det (l - (1 - XE^n-j.nXEY 1 XE^n-j ® Ipn-jXE) 

= det (1 - XE^n-j,nXE) ■ [l ~ ((1 ~ XE^-i.nXE) -1 $n-j, 1pn-j)E] (75) 

= det (1 - XEK n -j,nXE) ■ [1 - (^n-j,^n-j)E 

~ ((1 - XeKu-j^Xe)' 1 XEKn-j^XE^n-j^n-^E] ■ 

The advantage of using this formula is that VVi-j appears in the inner product (?/>n-j, VVi-j) e 
and the function (K n ^j ;n XEipn-j)(%)- We will see that we have good estimates on both of these 
quantities due to the fast decay of ip n -j(y; n) and K n ^j tTl (x, y) as n — > oo for all y £ E. 
We study two kinds of intervals E = 1% and E = j£(x*), x* > e. 

(a) For E = J„(a;*) where x*, which may depend on n, is in a compact subset of (e, oo), 
from (346), det (l — XjT( x *)K n ^j >n XjTr x \) — > 1. For the last inner product in (75), note 
that the operator norm of (1 — XjT^x^^n-j^XjT^,)) -1 1S uniformly bounded from (344), 
and 4) n -j{x) — > in L 2 (J^(x^,)) from (333). We will show that the -L 2 (jJ(x*)) norm of 
K n -j t nXjT (x^n-j is uniformly bounded. Hence we will have 

F n -j +1>n (a; Jn{x*)) = 1 - (tjj n -j,^n-j)jT {x ^ + o(l) (76) 

Therefore, we need 

(i) asymptotic evaluation of (ip n -j,ip n -j) jti x \, 

(ii) uniform boundedness of L 2 ( jj (re*)) -norm of K n -^ n XjT( Xt )i ) n~j- 

(b) For E 1 = /J, from (345), det (l — Xi^^-n-j,nXiT) — > Fq(T). For the cases in Section [3J we 
will show that K n -j >n XiT'ipn-j is uniformly bounded in L 2 (/J) but for the cases in Section [5J 
we will see that the L 2 (/J)-norm of K n ^j^ n Xjrip n -j is (^(n 1 / 6 ). We here state the necessary 
estimates separately. 

(bl) For Sections [5J we need 

(i) asymptotic evaluation of {ip n -j, i>n-j)iT, 

(ii) uniform L 2 (lJ) boundedness of K n _j jn XiTip n -j. 

Then it follows that, since the operator norm of (1 — Xi T K n -j,nXi T )~ l is uniformly 



bounded from Corollary 6.3 (c) and ipn-j — >■ in L 2 (/J) from (332), that 



n-j 



+1 ,„(a; II) = (F (T) + o(l)) • [l - + o(l)) . (77) 



(b2) For Section g we need 

(i) asymptotic evaluation of (i/j n -j,ip n ^j)jT, 
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(78) 



(ii) asymptotics evaluation of 

u j,n(0 ■ = 
in L 2 ([T,oo)). 

Then since (1 - X[T,oo)^n-i,nX[T,oo)) _1 ~+ (1 - X[T,oo)^Air y X[T,oo)) _1 in operator norm 
from (|343| > and rrV^ n _ i ( e +/r 1 rr 2 / 3 £) ~%n( e ) Ai (0 ^ in L 2 ([T,oo)) by Corollary 



6.1 



d)j we find that if Uj )T 



in L\[T, oo)), then 



P„_ j+lin (a;/J) =(F (T) + o(l)) • [l - <^ n _„ V„-,)/t 

1 n (79) 

- -g#J>(e)((l - X[T,oo)^AiryX[T,oo)) «n> Al)[ Tj00 ) + o(l)J . 

In Sections [3| [4] and [5] we only consider a > 0. When a < 0, the largest eigenvalue in the 



spiked source model defined by (21 ) has the same distribution as the negative value of the smallest 



eigenvalue of the spiked source model that is defined by the same formula but with the potential 
function V(x) = V(—x) and the external source matrix — A n _j+i. Since V{x) is regular as long as 
V(x) is, and the non-zero eigenvalue of A n _j + i is positive, the analysis in this paper applies for 
that spiked source model. We need to keep track of the smallest eigenvalue in the new spiked source 
model, and it can be done in the same way that we analyze the largest one. It can be checked that 
the limiting distribution of the smallest eigenvalue is not affected by the positive external source 
eigenvalue o, corresponding to the a < case of Theorem 1.2;[a) We skip any further remarks. 



3 When < a < \V{e) 

As outlined in Section [ij we need to show that the 1? norm of K n -j^ n xiTtp n -j is uniformly bounded 
in n, and need to evaluate (VVi— j) ipn-j)E asymptotically for appropriate choices of the interval E. 
The key part is the asymptotic evaluation of the function ip n _j(x). In Subsection 
evaluate T n _i(a) asymptotically and then use this in Subsection 



3.2 



to evaluate tp n - 

remaining subsections are devoted to the proof of the main theorems in each sub-case. 



3.1 



we first 
j(x). The 



3.1 Asymptotic evaluation of T 



n-3 ^ 



n-j 



[a;n) 



From the definition (70), 



/c poo 
-oo J c 



(80) 



for any cel. We take c = c(a) as in Definition 1.1 Note that c(a) > e since a < iV'(e). The 



reason that we split the integral at y = c will be clear in the below, particularly the paragraph 



before (|89j). 
Let 



(Cf)(z) :- 



1 



2ttz J R y- z 



dy 



(81) 
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denote the Cauchy transform of function / G L 2 (M). Using C+ — C_ = 1 and noting that <p n —j is 
analytic, we have 



^ 3 {y)e na Uy 



(82) 



Note that {C(p n -j){z)e naz — > exponentially as — > — oo since a > 0. Therefore, we can deform 
the contour and obtain 



where, with a constant Cr > 



{C^ 3 )(z)e naz dz, 



r+ur_ 



T + :={c + it | < i < C r } U {c + iC r - t \ t > 0}, 
r_ := complex conjugation of r + . 

The contours are oriented as indicated in Figure [5j Therefore we find 

r n _» = -/ (C< Pn „ j )(z)e naz dz+ ifn^e^dy. 
Jr + ur_ Jc 



(83) 



(84) 



(85) 



e = a N+1 



c(a) + iCr 



c(a) 



da) — iCr 



Figure 5: The contours r + and T_ 



Let 5 be given in Proposition 6.1 Let e < min{c(a) — e, 25} be a small enough positive constant, 
independent of n, such that all maximizers of G(x; a) in [c, oo) are in the interval (e + e, oo). Recall 
the asymptotics of {C(p n -j){z) summarized in Section |6j Since the contours T± lie in B$ (in 
Figure 10) in Section[6j from the asymptotic formula (320) for (Ctp n -j)(z), 



r+ur_ 



(Cifn-j^e^dz 



r+ur_ 



(86) 



Here we recalled the definition (32), H(z; a) := —g(z) + az + £. 
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We now use the method of steepest-descent to evaluate the integral asymptotically. By Lemma[TTTJ 
H'(c(a); a) = and H"(c(a); a) > 0. It is straightforward to check, with the help of the formula of 
g(x) in (14), that for z{t) = c + it, t > 0, the function ?RH(z(t); a) in t satisfies 

d 



dt 



MH(z(t);a) 



t 



Also for z(t) = c + iCr -t G T+, t>0, 

d 



dt 



&H(z(t);a) 



(c-s) 2 + t 



t — c + s 



{t-c + s) 2 + Cl 



V(s)ds < 0. 



$>(s)ds - a, 



(87) 



(88) 



is negative for all t > if C-p > l/(2a). Hence 5RH(z;a) decreases as z moves along T + counter- 



clockwise. Similarly KH(z; a) increases as z moves along T_ counterclockwise. Therefore T + U T_ 
is a curve of steep-descent for H with the saddle point at z = c. The fact that z = c is a saddle 
point of H is the reason that we have split the integral in (80) at c. 

From Proposition 6.1 [b) , Mj >n = ■Mj,n(z)(l + 0(n^ 1 )) and M.j >n (z) is analytic in B$. Moreover, 
M j>n , M' j n and 1/M jt n are uniformly O(l) in a neighborhood of c(a) and Mj )Tl — 0(z? 1 ) uniformly 
in n as z — >• oo. Thus the method of steepest-descent can be applied to (86) and we obtain 



f tr, no*; iV^Mj n ( c ) e »(H(c(«);a)-</2) 

/ (C^ n _,)(z)e™ a ^z = ^ (1 + o(l)). 



(89) 



Now consider the second integral in (85). From the asymptotics (318) for ip n -j in B$, 



^{y)e na Uy = \ M hn {y)e< G ^~^ dy 



(90) 



Using that Mj >n is uniformly bounded in any compact subset of [e + e, oo), and Mj tn (y) = 0(y 3 ) 
as y — > oo, and using that G(y; a) —> — oo at least linearly, we obtain the trivial estimate that 



(91) 



where G max (a) := max{G(y;a) | y > c(a)}. Together with (89), we obtain the following result. 



Recall the properties of G and H in Subsection 1.4 



than (91) and we obtain 



Suppose that a < a c . Then H(c; a) > G max (a). Therefore, (89) is exponentially larger 

(92) 



n-j 



e -nl/2 -iMj, n {c(a)) c nH(c(a);a) m , q(1)) 
y/H"(c(a);a) ' ' ' ' 



Suppose that a c < a < gV'(e) (assuming that V is such that a c < |V'(e)). Then G r 



a > 



H(c;a) and hence ( |9l[ ) is exponentially larger than (89). Suppose that a ^ Jy and let 
= %o(a) G (c(a), oo) be the unique point G max (a) is attained. If G"(xo;a) ^ 0, using the 
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Laplace's method applied to ( |90| ) (using the properties of Mj )Tt in Proposition 6.1 (a)), we 
obtain 



2tt 



n-j\ 



-nG"(xo; a 



-M j , n (xo)e nG( - xo ' a ^ n ^ 2 (l + o(l)). 



(93) 



If G"(x ; a) = and G( 2/c )(x ; a) ^ 0, G^')(x ; a) = for j = 1, • • • , 2fc-l, then the Laplace's 
method implies that 



r n -j(o) 



(2fe)! 
-nG( 2fc )(xo; a 



- x2fe dxA^ iin (xo)e nG ( xo;a )- n ^ 2 (l + o(l)). (94) 



When a € Jy, the contributions to (91) at each maximizer should be added. Examples of 
this case are in (155) and (174). 

If a = a c (note that since we assumed a < ^V'(e), this implies that V is such that a c < 
|y'(e)), then H(c;a c ) = G max (a c ). If we further assume that a c ^ Jy and G"(xo(a c ); a c ) / 
0, ((891) and (91) are of same order. Then by Laplace's method applied to (91), 



_,(a c ) = J-e-^ • A ^" ( "° (ac)) + " ^" (c(ac) ) e «G(*o(-)^)(l + (1)). (95) 
J V n l v ^-G"(x (a c )) V^W) / 



We can consider a double scaling case when 



a = a c + 



(96) 



where a is in a compact subset of R. By the definition of c and xo, a direct computation 
shows that 



G(x (a);a) =G(x (a c ); a c ) + ^ + 0(n" 2 ), 



H(c(a); a) =H(c(a c ); a c ) + 



aca. 



+ 0(n~ 



(97) 
(98) 



This implies that 



r„_j(a) = 

^ n£/2 ( M j:n (x (& c )) 



+ 



-iMj >n (c(a c )) 



V-G"(x (a c );a c ) V H "( c ( a c); a < 



; a(c(a c )-xo(a c )) j e «G(i (a);») /j _|_ (l)), 



(99) 



If G"(xo(a c )) = 0, then the term (91 ) is greater than (89) by a fractional power of n (see (94)), 
and hence the term involving A / fj i?1 (c(a c )) disappears in the expression of r n _j (ac). On 
the other hand, if a c 6 Jy, then there are more than one maximizers of G(y;a c ) making 
contributions in (91). We do not state the formulas explicitly here but instead state them in 
the appropriate subsections where they arise. 
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3.2 Asymptotic evaluation of ip n -j{x) := y) n _j(x;a) 
3.2.1 Algebraic formula 



From (71 ), 



J poo 

i=o J -°° 



(100) 



This can be written in the following way. Let c be any constant such that c > e. We will take 
c = c(a) as in Definition in the subsequent sections for asymptotic analysis, but the following 
result holds for any c > e. 



Lemma 3.1. For < a < ~V'(e), we have, with T± given in (84), for x G M \ {c}, 
r n _,(a)^(x) = e^-^)/ 2 )l (Cj0o) (x) 

/>oo 

K n _ Jin (x,z)e^^- / K n _ J> (x,y)e"^- v '^/ 2 )dy 



+ 



r+ur_ 



(101) 



where 



w , x ._ 7fc-i j)k{x)(Cip k -i)(z) - ^ fc _i(x)(Cy7 fc )(z) 
7fc £ - z 



(102) 



and Kk n (x,y) is defined in (69). 



Proof. By the same calculation that leads to (85), (100) equals 



T n . j {a;n)^ n . j (x)=e n{ - ax - y ^^ 
n-j-l 



n-j-1 



+ V Vi(a?) f (C^i)(z)e naz dz- V ^(s) / 
i=o ^r+ur- i=o ' /c 



(103) 



We exchange the sum and the integral in both terms. The second sum can be simplified by using 
the Christoffel-Darboux formula and becomes the last integral in (101). To analyze the first sum, 
we first take x £ C\ (KU T + U r_). By using the definition of Cauchy operator and from the 
Christoffel-Darboux formula, we have 



n-j-l 



n-j-1 



_i_ 1 /" 1 \ \ 

V Tpi(x)(C(pi)(z) = 7T-- I Y] ^i{x)<Pi{w)dw 

i=0 JK i=0 



Jn-j-1 1 /" i>n-j(x)<Pn-j-l(w) - i>n-j-l(x)<Pn-j(w) 



Jn-j 



2iri 



(w — z)(x — w) 



dw. 
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Using the partial fraction formula and the definition of the Cauchy transformation again, this equals 



Jn-j-l 1 1 
7n-j x — z 2iri 



1p n -j(x)ip n -j-l(w) - lj) n -j-i{x)<pn-j{w) 
W — Z 

t/j n -j(x)(p n -j-l(w) - 1p n -j-i(x)ip n -j(w) 



dw 
dw 



Jn-j-l 
ln-j 



Jn-j-l 



W — X 

1p n -j(x)(C<p n -j-i)(z) - 1pn-j-l(x)(Cip n -j)(z) 



x — z 



tp n -j(x)(Cip n - j - 1 )(x) - ip n -j- 1 (x)(C(p n -j)(x) 



x — z 



1p n -j(x)(C(p n -j-i)(z) - lf) n -j-l(x)(C(p n -j)(z) 



ln-j 

1 e -nV(x)/2 

2iri x — z 



x — z 



(104) 



where the identity (326) is used in the last line. Hence the first sum on the right-hand-side of (103) 
satisfies 

n-j-l „ 

V ViOc) / (api)(z)e naz dz 
•/r+ur_ 

1p n -j(x)(C<p n -j-l)(z) ~ 1p n -j-l(x)(C(p n -j)(z) 



i=0 
Jn-j-l 

ln-j Jr+ur 

P ~nV(x)/2 



e naz dz 



(105) 



x — z 



+ 



2iri 



-dz. 



r+ur_ x — z 



Note that this was proven for x € C\ (EU T + U r_), but the identity holds for x € M \ {c} as 
well by analytic continuation. The last integral equals — 2irie nax for x £ (— oo,c) and equals for 
x £ (c, oo) by Cauchy integral formula. Therefore we obtain (101). □ 



3.2.2 For x > e + e: 



Take c = c(a) as in Definition |1.1| in the formula of Lemma 3.1 Fix e > small enough so that 
[c, oo) C [e + e, oo). 



Lemma 3.2. For x £ [e + e, oo), 

= e "( G W-HW)/ 2 1 -^L - - e " H ^ 1 (C|00) (s) + 0(v^(l + M)^')} 
as n-^ oo and j = 0(1). 



(106) 



Proof. Fix e' £ (0, c — e — e). Assume that x satisfies \x — c| > e'. Noting ax — V{x)/2 
(G(x) + H(x) — £)/2, by Lemma |3.1| we have 



4>n-j(x) = e 



n(G(x)-U{x))/2 



{^7ra enHWl(C - )(X) + ^ (X) }' 



(107) 
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where 



r n _j(a) 



r+ur_ 



e nmz)-G( X ))/2 Kn _^ y)e n(ay-V { y)/2) dy 



From (336), the second integral over (c, oo) in (108) is 

Of (1 + |x|)" J ' ^ e n(G(y)-e/2)^ + \y\)-3 d y\ . 

On the other hand, for the integral over T + U T_, ( |318 ) and (320) imply that 

K n _ J> (x, 2 )e 7l ° 2 



(Mf |Tt (a:)Mf + i iTt (z) - Mj + i iW (x)Mj in (z)) w (G(g)-H(a:))/2 n(H(z)-l/2) 



a; — z 



By Proposition 6.1 [a][bj 



(M j>n (x)M j+ltn {z) - M j+ x !n (x)M j>n {z)) 



O 



x — z 



(i + \x\)-3\i + \z\y 

\x — z\ 



(108) 



(109) 



(110) 



(111) 



in (108) is 



Thus using the fact that c is the saddle point of H(z), we have that the first integral over T + U T_ 

(112) 



o((l + \x\)-i [ e n(H(z)-e/2) (l + \z\y d \ =o((l + bir^-VV 11 ^ 
V Jr+ur_ \x-z\ ) V J 



On the other hand, for T n _j(a), we have from (85), (86), (89) and (90) that 

r n _» = O (n-VV 11 ^)) + O QH e< G ^ 2 \l + |y|)-^) . 



(113) 



Note that we also have a matching lower bound. Comparing the estimate (113) of T n _j[a) and 
two estimates (109) and (112), we find that Qj^ n (x) = 0((1 + |x|) _ - ? ) uniformly for x > e + e if 
\x — c| > e' for a positive constant e'. Note that in this case the error term in ( 106 ) does not contain 
n. 



Now let x satisfy \x — c\ < e'. In this case, we start with the formula (101) with a different 
choice of c. We replace c by c± n -1 / 2 and let (r + ) =t U (F-) 1 * 1 be a contour deformed from T + U T_ 
by a semicircle of radius n -1 / 2 to the right/left, respectively, as illustrated in Figures gandQ Here 
we take the + sign if x — c > and take the — sign if x — c < 0. Then 



tpn-j(x) = e 



n(G(x)-H(x))/2 



^f-^e" H Wl (c±n _ 1/2i00) (x) + Qf n (o 



(114) 



and 



r n -j(a) 



(r + )±u(r_)J 



e n(H(x)-G(x)) Kn _^ n(X; z ) e "« dz 



c=pn 



-1/2 



e "(H(x)-G(x)) ijr ^ n ( a . >I/ j e n(oi/-V(i/)/2) dl/ 



(115) 
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The second integral has the same estimation as that in (108). For the first integral, note that |x — 
z\ > n -1 / 2 . Using this, and by recalling the asymptotics of a Cauchy-type integral /(r )±u(r )± nhz^ 
0(e nH ^) for such x, we find that the first integral of (115) is O ((1 + |x|) _ -?e nH ( c )) instead of 
O ((1 + |x|)~ : 'n~ 1 / 2 e nH ( c )) in the case when \x — c\ > e'. Hence we obtain (106) for \x — c\ < e by 
noting that e nli ^ /(e ni T n ^) = 0(y/n) for \x - c\ < n" 1 / 2 from ( pA3| ). □ 



(T- 



(T- 



(T- 



c + n 



-1/2 



(T-)- 



Figure 6: The deformed (T H 
(T_)+ for x-cG [0,e'). 



U 



Figure 7: The deformed (T_| 
(T_)~ for x - c G (-e',0). 



U 



3.2.3 For x near e: 

Let T be a fixed constant and let e be a small positive constant such that < e < min{c — e, 8q} 
where 5q is the constant in Proposition 6.1 and its corollaries in Section [6j Define the interval 



E Ttt := il \ (e + e, oo) = [e + p- l n~ 2 / 3 T, e + e] . 
For a given x G £^6) define £ by the relation 

x := e + ( S -1 n" 2 / 3 £. 
Lemma 3.3. We have for all < a < V'(e)/2, 

^ J (x) = 0(n 1 / e e-l^ /2 ), 

uniformly in x G Er,e and in n. 



(116) 



;ii7) 



;ii8) 



Proof. We use the formula (101). From (339) and (330), the integral over (c, oo) in (101) is 



O n^e-m 



||t|3/2 



a n{G{y)-l/2), 



l + \y\)->dy 



;ii9) 
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On the other hand, substituting (331) and (320) into (102), the integrand of the first integral over 
r+ U T_ in (101) is 

K n _ j:n (x,z)e naz = 0{n l l & e-iW\l + \z\y e^^^) , (120) 

for all x E Et,e an d z E T± since \x — z\ > c — e — e > 0. Thus the the first integral over T + U T_ 
in ((Toil) is 

() ( n^e-l^ 372 



I e n ^' £ / 2 \l + \z\ydz) =0 
Jr+ur- J 



n l/6 e -l\^ n -l/2 e nH(c) 



121) 



Substituting (119) and (121) into (101) and noting that l( Cj00 )(ic) = for x G £r,ej we obtain 



r n _,(a)^_,(x) = nVSe-IISI 3 / 2 (n" 1 ^) +o( e < G ^ 2 \l + |y|)-^y 



Comparing with (113) as in the previous subsection, we obtain (118). 



122) 
□ 



3.3 Proof of Theorem |l.]||(a)| and Theorem |1.2j|(a) 



Recall the outline of the proof described in Section [2j The proof proceed exactly same for both 
convex and non-convex potentials. The only important assumption is that < a < a c . 

We first evaluate ip n -j(x). Fix < e < So to satisfy the conditions in Subsection 3.2 where So is 
the constant in Proposition |6.1| and its corollaries in Section [HJ Since < a < a c , the asymptotics 

1 



(92) implies that 



- e "H(x) =0 (^ e «(H(*)-H( C )h 



(123) 



Since H(x) > H(c) for all x > c and H(x) — > oo fast by Lemma |1.1[ this term is larger than 
0[y/n(l + | | ) J ) . Inserting this into ( |106 ), we obtain 



tpn-j{x) 



' O (0^n(G(*)+H(*)-2H(c))/2) > x > ^ 

O (y/ne"^^)- 11 ^))/ 2 ) , e + e < x < c. 



(124) 



On the other hand, for x G Et )£ := [e + - — T ^ 2/3 , e + e] (see (116)), we have from Lemma 



3.3 



^ 3 {x) = 0{n 1 / & e-'^\ 3 ' 2 ) 



that 
(125) 



where £ is defined by (117). 



Now evaluate the inner product (tpn-j, ' l Pn-j)i T - We divide the interval 1% into two parts: 
(e + e, oo) and Er,e- From the asymptotics (124) of if) n -j and (330) of ip n -j, 

(^- i ^n-i)(e+e 1 oo)= f 0(^e n ( G M- H ^ (1 + 
Je+e 



+ y 0(v^e n(G(:c) - H(c)) (l + |x|)- J ')dx 
=0(e- e ' n ) 



(126) 
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for a constant e' > since when a < a c , G(x) < H(c) for all x > c and G(x) — H(x) < for all 
x > e (see Lemmas 1.1 and 1.2). On the other hand, by the asymptotics (125) of i/j n -j and (331) 



of ipn-j, we find, after the change of variables x t- > £ defined in (117), that 



T 



pn z l A 



(127) 



Therefore, 



(^_ J ,Vn- i >/T=0(n- 1 /3). (128) 

Finally, we show the uniform boundedness of K n ^j^ n XjTip n -j in L 2 (I^). From the asymptotics 
of K n ^ given in Corollary p.2| and the asymptotics ( |124[ ) and ( |125[ ) of ip n —j > we find for x E -Er,e 
that 



n-j 



)(*) = / 

Je 



2e/3n 2 / 3 



K n _ jin (x,y)il) n -j(y)dy+ I K n _ jtU (x,y)ij n - j (y)dy 



e+2e 



0{n 2/3 e -l\^-l\^\i/e e -l\v\^ ) ^ 
+ f 0(n 1 /6 e -fl«l 3/2 e «(Gfo)-H( J/ )) (1 + | y | rJ) ^ 

Je+2e 

/oo 
0(n l/6 e -||«| 3 / 2 e n(G( y )-H( C )) (1 + lyl) -J )dy 

-0{n^e-l^ 12 ) 



(129) 



where £ is defined by (117). Similarly, for x > e + e, 

(K n - jin x I T'ip n -j)(x) = / K n - j>n (x,y)4> n -j(y)dy + / K n _ jyn {x,y)ip n - j {y)dy 



T,e/2 
§e/3n 2 / 3 



e+e/2 



+ 



e+e/2 
oo 



( n l/6 e n(G(x)-H(x))/2 e -f|^|>/a n i/6 e -||,|»/aj_rf»7_ 

/3n 2 / 3 

0(v^e n ( G (*> -H (*^e n ( G M- H M))dy 



(130) 



0(^e"( G ^- H W)/ 2 e m ( G ^- H ( c »)d 2/ 

:O / e n(G(0!)-H(x))/2\ 



From these, we find that 

1 1 K n-j,nXlT 4>n-j \ \ \t. (jT) 



Hence from (77), we obtain 

lim det (1 - Xi? K n -j+i,nXiT 
Theorem l.]pij and Theorem 1.5 a) are proved. 



(131) 
(132) 
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3.4 Proof of Theorem 



when a < ^V'(e) 



The proof of Theorem l.^£b) is divided into three cases, a < |V r/ (e), a > ^^'(e) and a = ^V'(e). 
The first case is in this subsection. The second case is in Section HI The third case is discussed at 
the beginning of Section [5j 

We assume that a G (a c , hV (e)) and a £ Jy. Let xq = xq(o) be the unique maximizer of 



G(x) in (c, oo) as in Lemma 1.3 We assume that G"(xq) 7^ 0. See Remark 3.1 at the end of this 



subsection for a discussion when G"(xo) = (see (|59|)). 

Recall the definition of the interval Jj(xo) in (24). Note that Jj(xo) C (c, 00). We first evaluate 
tp n -j(x). By using (93), we have 



1 n-1 



-e nH(l) = O(V^e n( - n{x) - G{x0)) ), x>c. 



(133) 



Lemma 1-lptJ and |(b)| imply that H(x) increases monotonically in x > c and H(x) > G(x) for 
all x > e. Hence there exists e' > such that H(x) > G(xo) for all x > xq — e'. In particular, 
H(x) > G(xo) for x E J^{xq). Therefore (106) yields, noting that H(x) 
Lemma l.l^e) 



00 fast enough by 



4>n-j(x) 



B n(6(*)-H(*))/2. 



(l + o(l)), xEjJ(x ). 



(134) 



n-j\ 



Inserting the explicit asymptotics (93) for T n _j(a) into (134), we have for x > xq — e' where e' is 
the positive constant mentioned above, and in particular for x E jj(xo) that 



n-j\ 



-nG"(x ) 



1 



2vr M jt n(x ) 



,n(G(aO+H(a!)-2G(xo))/2, 



l+o(l)). 



(135) 



The inner product (ip n -j,^ n -j)jT^ X0 ^ is evaluated by using (135) and (330). For x E J„ (xo 

i) n - j {x)tp n . j {x) = 



-nG"(x ) M^ n (x) d „ (g(i) - G(io)) ^ + o(1 ^ 



2vr 



(136) 



•Mj, n (a;o) 

From the assumptions for the Theorem |1.2[ G(x) in (c, 00) has the unique maximum at x = xq 
and G(z) = G(xo) + |G"(xo)(x - xo) 2 + 0(|x - x | 3 ) for x close to xo where G"(xo) < 0. Also 
Mj >n (x), M'j n (x) and 1/Mj >n (x) are bounded uniformly in n for x in a compact subset of (e, 00) 
and Mj >n (x) = Aij jTl (x)(l + o(l)) from Proposition 6.1 [a),. Hence the standard Laplace's method 
applies and we obtain 



(i>n-j,ipn-j)jT( X0 ) = -7== / e-^ 2 d£,+o{l). 



(137) 



We now show that K n ^j in XjTr Xo \Tp n -j is uniformly bounded in L 2 (jJ(xo)). From (135) and 
the part (a) of Corollary 6.2, for x E JJ(xq), 



(K n - j} nXjT ixo) $ n -j)(x) 

-- [ 0(V^e n(G ^- H(x »/ 2 (l + \ x \)-j e n(G(y)-G(x ))^ + 
Jj ? T(x ) 

:0(e n(G (a; )-H W )/2 (1 + ^-jy 



(138) 
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Therefore, 



for some e' > 0. 



\ K n-j,nXjT( X0 )1pn-j\\L2(jT( X0 )) = 0(e '"), 



(139) 



Therefore, from (76), we obtain 



Jim det (l - X jT {xo)) Kn- j+ i,nXjT {xo)) ) = 1 - -j= J e~^d^ = G(T), 



(140) 



and Theorem 1.1 [b) is proved. 



Remark 3.1. When G"(xq) = 0, the Gaussian function e 2^ in (137) is replaced by a higher-order 



function such as e ^ 2k (k > 1). The rest of the proof is very similar. The result is the limit theorem 



as in (59). 



3.5 Proof of Theorem |1.3j|(b) 



Let V be a potential such that a c < ^V'(e) and a c Jy. We assume that G"(xo(a c ); a c ) / 0. Let 



a = a c + 



Q: 



141) 



n 



where a is in a compact subset of R. 

First, consider ip n _j(x). Note that the estimates (124) and (125) still hold. However, G(xo(a c ); a c )+ 
H(x (a c );a c ) - 2H(c(a c );a c ) = at x = xq, so when a = a c , (i> n -j, if> n -j) (e+e,oo) is 110 longer ex- 



ponentially small. We need an asymptotic formula of ip(x) like (135). By inserting the asymptotics 



(99) of r„_i(a) into (106), similar to (135) we obtain for x > xq — e' where e' is a positive constant 



defined similarly as the e' in (135), and in particular x G J^, that 
Ci(a) /-nG"(i (o)) 1 



Ipn-j(x) 



n(G(x;a)+U(x;a)-2G(x {a);a))/2 , 



l + o(l)), (142) 



C + Ci(a)V 2vr Mj, n (x (a)) 
where (we omit the dependence of Co and C\(a) on n and j to make the notations simple) 

-iMj in (c(a c )) 



Cn 



V / H"(c(a c );a c ) : 



c / x = jWwWa^ e ^(a c )-c(a e )) 
V / -G"(x (a c );a c ) 



(143) 
(144) 



The constants Co and Ci(a) are positive from Proposition 6.1 If a is fixed, Co, Ci(a), C 1 and 
Cf x (a) are uniformly bounded in n. Set 



Co 



C + C!(a) 



(145) 



From the definition, pj^ n (a) is a decreasing function in a, pj ;n (a) — > as a — > 00 and p Jin (a) — > 1 
as a — > —00 for each fixed n. Also for a fixed a, pj^ n (a) is in a compact subset of (0, 1) uniformly 
in n. Note that when the support of the equilibrium consists of one interval, then Cq and C\{a) 
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are independent of n, and hence so is pj^ n {a). We prove formulas (48) and (47) in Theorem 1.3 
separately. 



The proof of ( 48 ) is similar to that in Subsubsection 3.4 For the evaluation of (ip n -j > ipn-j) jj(x ) > 
we repeat the arguments of (136)-(137). Noting that the formula of ip n -j in (142) is the same as 



that in (135) except for the multiplicative factor C\{a) / {Cq + C\(a)), we obtain similar to (137) 
that 

/OO 
e~5? 2 -Ji=(l + (l)), (146) 

,(°)/ 



3.4 



where p)n( a ) is defined in (145). The estimate (139) follows from the same calculations in Subsec- 
tion 



and we obtain from (76) that 



det (l - XjT( X0 ))Kn-j+l,nXjT( X0 ))) = Pj,n( a ) + i 1 ~ PjA a )) G ( T ) + 



(147) 



and (48) is proved. 



We now prove (47). When a is given by (141), the estimate (124) still holds. Similar to (146), 



we obtain by estimates (142), (124) and (118) that 



/OO 
-=e-^(l + o(l)) = (1 - Pj> (a))(l + o(l)). (148) 
-oo V27T 

To show that K n -j iTl XiTip n -j is uniformly bounded in L 2 (/J), we proceed as in (129) and (130), 
and obtain 



{K n -j,nXlT^n-j){x) 



OinWes 151 7 ), x E E T>e 

G ( e n(G(x)-H(z))/2) x > e + e 



(149) 



where £ is defined by (117). Hence \\K n -j jn XiT'4>n-j\\L 2 (I T ) = 0(n 1 / 6 ). Therefore, by (77), we 
obtain 

det (l - XiTKn- j+ i,nXiT) = Pj , n (a)F (T) + o(l) (150) 



and (47) is proved. 



3.6 Proof of Theorem 



1.4 



when a E (a c , \V{e)) 



We prove Theorem 1.4 when ao < iV(e). The case when ao > \V{e) will be discussed in 



Sections |4] and [H 

Let ao E (a c , ^'(e)) and ao E Jy . Hence ao is a secondary critical point. In this case, the 
maximum of G(x; ao), x E (c, oo), is attained at more than one point. The case when the maximum 
of G(x; ao) is attained at more than two points can be attained by a straightforward extension and 



this yields (52). We omit the details in that case. 



Denote the two maximizers of G(x;ao) by x\ := x\(ao) and x<i := ^(ao)- Let xi(ao) < ^(ao)- 
Assume that 

G"(a;i(ao); oq) + 0, G"(x 2 {a ); a ) + 0. (151) 
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The case when one of the derivative vanishes is discussed in Subsection 13.71 Let 

a 



a = ao + 



n 



(152) 



where a is in a compact subset of R. 

First we evaluate ip n -j(x). The asymptotics of T n ^j(a) in this case is not explicitly computed 
in Subsection |3.1| hence we first compute this by extending the formula ( 93 ) (see the paragraph 
following (94)). There are two differences from the case leading to (93) in Subsection 3.1 The 
first is that there are two maximizers of G(x;a>o) and the second is that a scales in n as in (152). 
The first difference simply results in adding the contributions from the both maximizers since both 
term are of the same order due to the condition (151). Regarding the second difference, note that 
since G(x;ao) has maximum at xi(ao) and 2:2(00), G(x;a) has two local maxima at two points, 
denoted by x\(a) and 2:2(0), which are close to xi(ao) and 2:2(00), respectively. (Indeed, one can 
easily check that X{(a) = Xj(ao) + _g"( "(a ))n 0{n~ 2 ).) Using the definition of G and the fact 
that Xj(ao) is a critical value of G(x; ao), we find 



d_ 

da 



G(xi(a);a) = Xi(ao) 



1,2. 



a=ao 



Hence 



a 



G(xi(a);a) = G(xi(a );a ) + Xi(a )- +0(n ), 

n 



1,2. 



(153) 
(154) 



Therefore, as in (93) we obtain as n — > 00 (note that G(xi;ao) = G(x2;ao) and xi := Xi(ao)) 

-M hn { Xl )e^ a + 



2tt 



2tt 



-nG"( Xl y 



-nG"(x 2 



-M ])n {x 2 )e x * a 



1 + 0(1)). 

(155) 



With this asymptotics of T n _j(a), the rest of the analysis is similar to (135), and we obtain for 
x G Jn(xi), 



2tt 



-nG"(xi) 



Mj, n {xi)e x ^ + 



2tt 



-Mj, n (x2)e 



-nG"(x 2 )' 

x e n(G(a;;a)+H(a;;a)-2G(a;i;ao))/2/j + n\\ (156) 

We now compute the inner products (tjj n -j,tpn-j) JJ^,)' * = 1?2. Using (330), (cf. ( |136 )) 



i> n -j(x)ifj n -j(x) 



■ WnjlJe -(1 + 0(1)) 



A 1 (a) + A 2 (a) 



where 



Ma) :-- 



M jin (xi(a );ao) pXt{ao)a 
y / -G"(x i (ao);a Q ) 



1,2. 



(157) 



(158) 



Like Cq in (143) and C\(a) in (144), Ai(a) is positive and is of finite distance away from uniformly 

.'■ - - . j :. (159) 



in n. For each i = 1,2, if we set 



y/-nG"( Xi y 
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then for £ in a compact subset of R, we find using the Taylor expansion in x, and (154) that 



G(a:; a) = Gfa; a ) + x t - - |- + 0(n 3/2 ) 

n in 



Thus we find 



ip n -j(x)ip n -j(x) 



nG"(xi;a ) Aj(a) 
2-7T Ai(a) + A 2 (a 



■e-^ 2 (l + o(l)) 



(160) 



161) 



for x given in (159) and £ in a compact subset of R, for each i = 1,2. Together with an easy 
estimate when x is away from x±(ao) and £2(00), this implies, as in Subsection 3.4, that 

^n-j,4>n-j)jT {x2) =p ( ^ n {a) 
^n-j^n-j)jT {xi) =pfl{a) 

for any fixed T, where for i = 1, 2 





|€ 2 












r 


7t V2tt 




Ai (a) 





+ o(l) 



The properties of p^(a) stated in Theorem 



A 1 (a) + A 2 (a)' 

can be easily checked. 



(162) 
(163) 

(164) 



1.4 



The L 2 norm 1 1 K n _^ n XjT ( Xn )tp n -j \ \ L 2 / jt ( x .)) is estimated by the same argument as that for (139) 
above in Subsection . 



3.4 



by using the asymptotics ( 156 ) of ip n -j- The result is the same exponentially 



decaying bound. 



Therefore, by (76), we obtain 



det (l - XjT( Xl )K n -j+i,nXjT( Xl )) =p ( ll(a)G{T) + o(l), 

det (l - XjT( X2 )Kn-j+i,nXjT( X2) ) =pfl{a) + pfl(a))G(T) + o(l). 



(165) 
(166) 



Thus Theorem 



1.4 



3.7 Proof of Theorem 



when a < 2 U'(e) is proved. 

when a G (a c , |V'(e)) 



1.6 



We prove Theorem 1.6 when ao < ^V'(e). The case when ao > ^V'(e) will be discussed in 
Sections H] and [H 

Under the assumption of Theorem |1.6[ for some k > 1 



G"(a:i(ao);ao) 

G (i) (x 2 (a );a ) =0 for % = 1, . 
G( 2fc )(x 2 (a );a ) ^0. 

We consider the double-scaling situation when 

logn a 



,2k- 1, 



(167) 



a = ao 



+ 



n 



where g := 



1 i_ 

2 2fc 



x 2 (a ) - xiiflo)' 



(168) 
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for a in a compact subset of M. 



The analysis is similar to Subsection 3.6. For each % = 1,2, we have, as in (154), 

-q log n + a 



G(xj(a);a) = G(xj(a );a ) + Xj(a )- 



n 



+ o(n 1 ). 



For 



x =xi(a ) + 



6 



yJ-nG'^x^aQ^aoY 



we obtain, as in (160) 



\ „, , s , , . — <?logn + a £? , _i N 

G(x;a) = G(xi (a );a ) + £i(a ) — + o(n ) 

n zn 



for £i in a compact subset of R. Similarly, for 



x =x 2 (a ) + 



(2*)! 



-nG( 2fc )(x 2 (a );a ) 



l/(2fc) 



6, 



we have (using (167)) 



G(x;a) = G(x 2 (a );a ) + x 2 (a )- 



-glogn + a £f fe _ x 
+ o(n ) 



for £ 2 in a compact subset of M. Therefore, as in (155) above (cf. (94)), 



2tt 



-nG"(xi) 



+ 



(2k)\ 



\ 1/(2*) 



G( 2fc )(x 2 ;a ). 



-n 



-00 



(l + < 



Since 



1 



1 



X2 XI 

+xiq = 7 - + x 2 q = J ^ 

z zfc x 2 — x\ 



we can write (174) as 

e ne / 2 T n ^(a) = n -^T e " G ^«) [ Bl (a) + B 2 {a)] (1 + o(l)) 



a=2 _ 51 
"2" 2fc 



where 



Jb^ (xi)eX1Q ' 



(2A;)! 
-G( 2fc )(x 2 



\ 1/(2*) ^oo 

— M^ n (x 2 )e x ^ e-^^ 2 ,. 
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As in (156), for x E JT(x\) 



~T 2k e 



"-"fir p r l (G(x;a)+H(x;a)-2G( a; i(ao);ao))/2 



B l (a) + B 2 (a) 



From this we obtain, as in (161), that 



^(^(x) =pg(a) A / - raG l Xi;a0) e-^(l + (1)) 



2vr 



for x given in (170) and £i in a compact subset of R, and 

.(2) 



ifj n - j (x)ifj n - j (x) =p\J i (a) 



-nG (2k \x 2 ;ao) 
(2k)) 



l/(2fe) 



e *2 



for x given in (172) and £i in a compact subset of K, where for i = 1,2 

(i) / n _ ^(") 
^' nlaj - Si(a) + fi 2 (a)" 



From the definition, the properties of p^p n (a) in Theorem 1.6 follow easily. 
Thus it follows as in ( 162 ) and ( 163 ) that 



1 /"OO 2£ 

$n-j(x),^n-j(x))jT {x2 . h) =pfl{a) / e-t 2k d£ + o(l), 

J-oo e ^ T 



(ipn-j(x),i'n-j(x))jT {xi) =p^l(a) J -^=e a^df + pJJ(a) + o(l). 



OO 

OO J 



.(2)/ 



(179) 



(180) 



(181) 



(182) 



(183) 
(184) 



As in the proof of Theorem l.^[b) in Subsections 3.4 and in the proof of Theorem 1.4 



Subsection 



3.6 



m 



we have ^Kn-^nXjT^.^n-j^jT^ -> 0. 



Thus we obtain, from (76), that 



det (1 - XjT( Xl) K n _ j+hn XjT {xi) J =pfl(a)G(T) + o(l), 



det (l - XjT {x2 . k) Kn- j+ i,nXjT (x2 . k) ) =pfl(a) + vfl(ot)) 



i-oo eS d $ 

/OO (:2k j/- 
-oo e? 



(185) 

+ o(l), (186) 



and Theorem 



1.6 



when ao < \V'(e) is proven. 



4 When a > |V(e) 



Note that if a > ^^'(e), then a > a c . In this section, we prove Theorem 1.1 b) 

It// 

2 



rems 



1.2 


», 


1.4 


and 


1.6 



and Theo- 



for the case when a (or ao) > 2^'( e )- After a small change at the first 
step, the analysis is the same as in the case when a c < a < ^V'(e) discussed in Subsection 3.4|3.6 



and 3.7 The proof of Theorem l.l||(b) is identical to the proof of Theorem 1.2 (b) 
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Note that c(a) = e in this case (see Definition 1.1). Since G'(e) > when a > ^V'(e), 
G max (a) := max{G(a;;a) : x G [e, oo)} satisfies G max (a) > G(e) = H(e). Let e > be small 
enough so that all the maximizers of G are in (e + 2e, oo) and 



G max (a) > H(e + 2e) > H(e + e). 
We have the following formula of ip n -j(x). This is the analogue of Lemma 
Lemma 4.1. Let a > iV'(e). As n — > oo while j = 0(1), 



(187) 



3.2 



n-j 



(a) 



ip n ^(y)e na ydy(l + o(l)) 



and 



e^Y n . 3 (a) 



3 nH ^ +0((l + \x\)- j ) 



(188) 



(189) 



for x G [e + 2e, oo) 



Proof. As in (85), we write r n _j(a) = T n ^j(a;n) as 



■ n-j 



(a) 



(Cvn- 3 ){z)e naz dz+ / <p n Jy)e na ydy, 



r + ur_ 



where, for a large enough but fixed positive constant Cj, (cf. the contour T defined in (84)) 
f + := {e + e + it \ < t < Cf } U {e + e + iCf - t \ t > 0} 



(190) 



191) 



and r_ is the reflected image of T + about the real axis. The contours are oriented as indicated in 
Figure [8l 



As in (86), by using (320) for (Cip n -j)(z), the contour integral over T in (190) satisfies 



r+ur_ 



(Ctpn-j^Z^dz 



r+ur_ 



Mj, n {z)e n ^ z ^ l ^dz. 



(192) 



Now H(z) has no saddle point in T± since H(z) ^ for all z £ T. However, it is easy to check 
that for z(t) = e + e + it, < t < Cf and z(t) = e + e + iCf — t, the formulas (87) and (88) still 
hold verbatim except that c becomes e + e, provided that Cf is large enough, say, Cf > ^. Hence 
5RH(z) decreases strictly as z travels along r + in the direction of the orientation. Similarly, 5?H(z) 
increases strictly as z travels along F_ in the direction of the orientation. Noting that for fixed n, 
Mj >n (z) — > zi~ l as z — > oo and for fixed z, Mj >n (z) is uniformly bounded in n from Proposition 
we obtain 



6.1 



/ (C^ j )(z)e naz dz = 0( e nH(e+e)-n£/2-j 

ir A ur_ 



On the other hand, consider the second integral in (190): 

/•oo 

Vn-AyY nay dy= \ M hn (y)e< G ^- l ^dy. 



(193) 



(194) 
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e + e + iCr 



-0- 



• e + e — iCi 



Figure 8: The contours T + and T_ 



By using the Laplace's method, we find an estimate similar to (91). Hence we find that (194) is 



exponentially larger than (193) due to the assumption (187). Thus (188) is proven 



Now consider i/j n _j(x). Analogous to ( 101 ) in Lemma 3.1 we have, for x 6 M. \ {e + e}, 



(a)4> n -j(x) 



n-3 

+ / K, 
Jr+ur_ 



K, 



n—j,n 



x,y)e niay - viy)/2) dy. 



(195) 



e+c 



Using this formula, due to the property of the H(x; a) on T± and G(x; a) on (e + e, oo), the analysis 
of the proof of Lemma |3.2| applies without any changes. If we restrict x > e + 2e, then the error 
term 0{y/n(l + |:c|)~- J ) in (106) can be replaced by 0((1 + (xl) - - 7 ) since \x — z\ > e for z £ T± as 
in the first part of the proof of Lemma |3.1| We skip the details. □ 

Since tp n -j is the only term that depends on a and its asymptotic formula for a > ^V(e) is 
same as the case when a c < a < ^V'(e) in Section [3J all the analysis in Section 3|hold without any 
changes. Therefore, we obtain the proof of Theorem |l.l]|^b)| and Theorems L^J)J 1.4 and 1.6 for 
the case when a (or ao) > |V(e). 



5 When a = ^'(e) 



First, suppose that a = ^V'(e) > a c . Then G max (a) := max{G(x;a) | x G [e, 00)} satisfies 
G max («) > H(e) = G(e) (recall Definition 34 and (33)). This property is enough to prove 



Lemma |4.1| and the analysis of Section [4] applies without any change. Hence we obtain the proof 

I V (e) > a c . Combining the results of the previous two 



of Theorem 



1.5 


>)■ 


1.4 


and 


1.6 



when a 



sections, we have proved all theorems except for Theorems 1.1 Kb) l-3j[a) and 1.5 
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Theorem l.l [b) and Theorem l.c a) share the same proof and this is given in Subsection 5.1 



The proof of Theorem 1.5 is in Subsection 5.2 



5.1 Proof of Theorem [TT](b) and |1.3|Ka) 



Let V be a potential such that a c = ^V'{e). We assume that a c ^ Jy. (This holds under the 
assumption of convexity of Theorem |l.l[ a).) Then G(e;a c ) > G(x;a c ) for all x > e. We consider 
a double-scaling situation when 



a c + 



f3a 



n 



1/3' 



(196) 



where a is a real number in a compact subset of ' 



5.1.1 Computation of r n _ 



Lemma 5.1. We have 



where 



r„_,-(a) 



Qn 



e° 3 / 3 (% ri (e) + (l)), 



Qn = Qn{a) := e 



n(G(e;a)-i/2) 



(197) 



(198) 



Here Bj n (z) is given in (322). Note that B, n (e) is in a compact subset of (0, oo) independent of n. 



Proof. In the proof of Lemma 3.1 when a < lV(e), we have taken the contour T± to pass the 



point c = c(a) at which 3tH(z;a), z £ r±, takes its maximum (see (|85|)). Near this point, we had 
H(z; a) — H(c(a);a) ~ «(z — c(a)) 2 for some constant k > 0. This quadratic term changes when 
a = a c = iV'(e). In this case, c(a c ) = e, and (note (30)) 



H'(z;a c ) 



-g'(z) +a c 



-g'(z) + g'(e) 



(199) 



vanishes at z = e. Now since g'(-z) is the Cauchy transform of the equilibrium measure fy(x) 
(see (11)), which vanishes like a square root at x = e, we find that there is a constant > such 
that H'(z; a c ) ~ n'(z— e) 1 / 2 for z near e such that z— e ^ R_. Hence H(z; a c )— H(e; a c ) ~ k(z— e) 3 / 2 
for a constant k > 0. This 3/2-order of vanishing implies that near z = e, 3?H(z;a c ) decreases 
most rapidly in the direction of angle 2-7r/3 and — 27r/3 as z travels away from e. 

With the above preliminary computation in mind, we defin^jthe contours E± as (see Figure |9| 

E + = E 2+ UE n US 3l S_ = SiUS 2 ,_US 2 -, (200) 



3 Here, the exact shape and the angle of the contour from z = e is not important. For example, we can use 
the contour that extends straightly upward from z = e as in Figure [5] with c replaced by e. The local behavior 
H(z;a c ) — H(e;a c ) ~ k(z — e) 3 ^ 2 near z = e shows that 5RH(z) decays as z travels vertically away from e at least 
locally. One can check 5RH(z) indeed decreases as z moves away from e along on the entire curve. Our choice of the 
contour E is made for the convenience of the formulas that appear later. 
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e + ujCy, 



CLN+1 



<e + ore 



^e + w 2 C s 



Figure 9: The contours S + and £_. 



and 



Si ={e + oj 2 e-iC s + t\ t< 0}, 
£ 2 ,_ ={e + oj 2 e + it | -C E < t < 0}, 
E 2 _ ={e - w 2 t | -e < i < 0}, 
£ 2+ ={e + wt | < t < e}, 
£ 2 , + ={e + oje + it\0<t< C s }, 

S 3 ={e + oje + iC s -t\t>0}, 



(201) 



where w := e 2?rl//3 . Here e is a fixed constant chosen to satisfy the condition (215) below, and Cs 



is a positive fixed constant large enough, say, greater than l/(2a). As in (85), we have 
r„-i(a) = -/ (Op n - J )(z)e naz dz+ I vn- 3 {y)e nay dy. 

Let 



s + us_ 



(202) 



£|± = S 2 ± n {z G C | \z - e| < n~ 11/21 }. (203) 

We first consider the part of the first integral in (202) over £|+- Inserting the asymptotics ( 324[ ) 
for (dp n -j)(z), there are two terms, one involving Ai(w 2 $(2:)) and the other involving Ai'(w 2 $(z)). 
We compute each of the integrals using the change of variables z t— )■ £ defined by 



2 := e + 



/?n 2 / 3 



This change of variables and the double scaling ( 196 ) imply that 



s n(-iV(*)+«) = Q ne ^(l + (1)) 



(204) 



(205) 
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uniformly in z E £§+, where (recall Lemma 1.1 ^c) ) 



g n _ e n(-iV(e)+oe) = g n(H(e;a)-£/2) = e n(G(e;a)-f/2) ^ ^06) 

as defined in (198). Therefore, using the property ( 310[ ) of $(^), and noting that \z — e| < n~ n / 21 , 
the two integrals involving Ai(uj 2 &(z)) and Ai'(w 2< l ) (2;)) satisfy 



^2+ 



Ai(o; 2 $(z)) J B iin (z)e n( -^ (2)+ ^ ) dz 



and 



/3n 2 / 3 



B i>Tl (e) Ai(£)e^ + (1)) (207) 



/ kl{uj 2 ${z))D j!n {z)e<-l v ^ + ^dz = |% (^, n (e) /°° Ai'(£)e^ + o(l) 
Jess /3™ 2/ci V 



(208) 



Observe that the integrals involving Ai(£) and Ai'(£) are convergent as these functions decay faster 
than exponential functions as £ — > +oo. From these, we find that 

Qn 



{Cy n -j){z)e naz dz 



B hn (e) J™ Ai(£)e<^d£ + o(l)j . (209) 



Now consider £2+ \ S|?P. By the property (310) of &(z), we have that for \z — e| < 1, there 
exists ci > such that 

|n~ 2/3 $(z) _ _ e )| < ci\z - e| 2 . (210) 

Hence the asymptotics of Ai(£) and Ai'(£) as £ — > 00 ([1, 10.4.59 and 10.4.61]) imply that z £ £2+) 
if |jb - e| < (3(1 - (3/4) 2 / 3 )/ci, then 



Hence (324) implies that 



Ai( W 2 d>(z))=0((l + |£|)- 1 / 4 e-^ 3/2 ), 
AiV$(z)) =0((1 + |£|) 1/4 e-5« 3/2 ). 

(C7^„_ i )Wet v W = 0(n 1 /6 e -^ /a ). 



Also for 1 2 — e| < 1, there exists c% > such that (cf. (205)) 



Que 



l p n{-\V (z)+az) 



0(e c 



jC+n~ 1//3 C2^ 2 ^ 



Hence if we take e in (201) small enough so that 

f3c^ (l - (3/4) 2 / 3 ) and e < /T 1 Q 



e < 1, e < 



-2 ' 



(211) 
(212) 

(213) 
(214) 

(215) 



then combining (213) and (214), we have, for n > (8a) /p 1 . 



L 

-Qn 
-Qn 



{Ccp n _ j ){z)e naz dz 
O f ri 1 / 6 e-^ 3/2 + a ^+ n ~ 1/3c2 « 2N ) 



(216) 
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For the rest of £+, by a direct calculation as in the inequalities (87) and 



we find that 



$l(H(z; a)) decreases strictly as z travels away from e+oje along X2'+U£3. Also by direct calculation 
we verify that 

e' := KH(e; a c ) - RH(e + coe; a c ) > 0, (217) 

where t is a positive constant depending on e. Since H(z; a) — H(z; a c ) = -^j^z — > as n — > oo for 
a fixed z, the difference (217) with a c replaced by a is also bounded below by \t for large enough 
n 



. Thus, from Proposition 6,l[b) 



L 



{C^ n . 3 ){z)e naz dz =0 



S 2 / , US3 



-l e n(mi(z;a)-e/2) d i. 



n(H(e;a)-£/2) p> I „- ine' - 



= e ™ ^0(e 
--Q n O(e-\ nl> ). 



(218) 



Combining p09| ), ( [216] ) and ( [218] ), we obtain 



(CV^X^e^cfe 



V" 



B iin (e) Ai(£)e^d£ + o{\)\ . (219) 



The integral over S_ can be evaluated in a similar way. Alternatively we can use the symmetry 
(C<p n -j)(z) = -(C(p n -j)(z). We have 



Q, 



Bj, n (e) Ai(0e^dt + o(l)) . (220) 



For the integral over (e, 00) in (202), we again consider three intervals (e, e + n~ n / 21 ], (e + 
n~ n / 21 ,e + e) and [e + e, 00), and proceed as before. We now use the asymptotics (322) for ip n ^j(z) 
in the first two intervals and ( |318 ) for the third one. Note the similarity of (322) and (324). The 
calculation is similar and we obtain 



Q, 



B j>n (e) Ai(0e a *d£ + o(l)\ . 



(221) 



Combining fl219|), (|220j> and (|221[), we find 

Q 



IV,(a) = ^= (bjM Ai(0(e Q « + e"* + + o(l)) 



But it is easy to check that 



J Ai(0 (e* + e<^ + 



d£ 



3 « 3 /3 



(222) 



(223) 



for a £ C. Indeed if we denote the left-hand-side of (223) by /(a), then by using Ai"(£) = £ Ai(£) 
and 1 + w + oj 2 = 0, we find f(a) = a 2 f{a). Recalling that J °° Ai(£)d£ = \ (see e.g. 0Q 9.10.11]), 
we obtain (223). Thus we obtain (197). 

□ 
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5.1.2 Evaluation of ip n ^j(x) 

We evaluate t[) n -j(x) for i 6 = [e I /3 _1 n~ 2//3 T, oo) as n — > oo. Let Jo be the constant in 
Section [U 



Lemma 5.2. Lei < e < 2<5o 6e the constant in (201), satisfying the condition (215). For 
x £ E Te/2 := [e + P~ 1 n~ 2 / 3 T,e + e/2], we have 



■{- 



C _ a ^) + e -« 3 /3( J_ e n(--V(x)/2) _ e< 



(224) 



where £ is defined by the relation x = e + /3 1 n 2 / 3 £ as in (117), Q n is given in (198) and C a (£) 

(225) 



is defined in (18). For x > e + e/2, we /icwe 



^ n _.( s ) = M^_^! e n(G(,;a) + H(,;a)-2G(e;a))/2 (1 +o(1)) _ 



Note that £ G [T, oo). Let C < T be a real number, and set C n := /3 x n 2 / 3 C. For x G Zj 

K n - j>n {x,z)e naz dz 



(hence x > e + C n ), we have, as in Lemma 3.1 
T n ^{a)^{x) = e »(«-n*)/2) + 



(S++C„)U(E_+C n ) 



K ri _, > (x,n) e <^- y ^/ 2 )dy. (226) 



e+C*„ 



Here S± + C n denotes the contour £± translated by C n . For example, £+ + C ra = (X24- + C n ) U 
(£2'+ + C n ) U (S3 + C n ), cf. (200). We divide the proof of Lemma 5.2 into two parts. 



Proof of (224). First we consider the integral over £+ + C n in (226). For x G £r,e/2 an d z £ 
£ 2 + + C„, from ([322} and ([3241), 



n 1/3 [/i(x, z) Ai($(x)) Ai(u; 2 $(z)) 

+ 14 (x, 2) Ai($(x))w 2 Ai'(o; 2 $(z)) 
+ Va(a?, z) Ai'($(x)) Ai(w 2 $(z)) 

+ n^l^x, 2) Ai'($(x))w 2 Ai'(w 2 $(z)) 



(227) 



e 7ri/3 e -nV( 2 )/2 



where 



Ui(x, 


z) = 


7n-j-l 


Bj !n (x)B j+ljn (z) 


- J B J+ i >n (x)£? J> (z) 






7n-i 


X 


— z 








z) = 


7n-i-l 


£? J> (x) J D j+ i jn (z) 


- B j+ljn (x)D j<n (z) _ 


Wt(x 




7n-i 


X 


— z 


x — 


z 


V2(X. 


z) = 


7n-i-l 


Dj, n (a;)Bj + i )n (2) 


- D j+ i jn (x)B jtn (z) _ 


W 2 (x 


,z) 


7n-i 


X 


— z 


x — 


z 


C/ 2 (x. 


z) = 


7n-i-l 


r^ n (x).Dj + i in (2:) 


- D j+ i }n (x)Dj >n (z) 







(228) 



7n- 



x — z 
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We now show that the main contribution to M. n -j >n (x, z) comes from the middle two terms on the 
right-hand side of p27| ). Clearly, U x (x,z) = 0(1) and U 2 (x,z) = 0(1). Since Wi{x,x) = 1 = 
— W 2 (x, x) from (329), it follows that 

Vi(x, z) Ai($(x))w 2 Ai{lo 2 <S>(z)) + V 2 (x, z) M'($(x)) Ai(w 2 $(z)) 
Ai($(x))w 2 Ai'(u 2 $(z)) - Ai'(*(x)) Ai(w 2 $(z)) 



x — z 



+ 



Wi(x, z) — W\(x, x 
x — z 

W 2 (x,z) - W 2 (x,x) , 



+ ^ 1 ' A\(Mx))u? Ai'(w 2 $(z)) 

x — z 



(229) 



x — z 



Ai'(*(z))Ai(u^*(«)). 



Observe that 



Wj(x,z)-Wi(x,x) 



= 0(1), i = l,2. For 



(230) 



using the estimates (211) and (212) for Ai(o; 3>(z)) and Ai (uj <i>(z)), and analogous estimates for 
Ai($(x)) and Ai'($(a;)), we find that (noting that the £ in (211) and (212) are slightly different 
from the £ and rj in ( 230 ) ) 



fin 



2/3 



^Ai^AiV^-Ai'tOAi^) + 1/ 3 e _i ( |^|3/ 2+|77| 3/ 2) 



>) 



e 7ri/3 e -nV( 2 )/2 



(231) 

for x and z in (230). For z S (S 2 '+ + C n ) U (£3+ + O n ), noting that \x — z\ > 1/e, a straightforward 
calculation using ( [322] ) and p20| ) implies that 

K n ^ n (x,z) = 0(nV*e-*W a/ \l + |^|)^ ft (-«W + «/ 2 )). (232) 

Hence we obtain, by noting (205), ( 214[ ) and the condition (215) satisfied by e, using —g(z) +1/2 + 
az = H(z; a) — ^/2 and noting the calculation in (218), that 



K n ^Jx,z)e naz dz 



-Qn 



si p 00 Ai(£) Ai^w 2 ^)^ 2 - Ai'(£) Ai(w 2 r?) 



e 3 



e ar >dri + 0{n-^ s e 



(233) 



The estimates for the integral over the contour £_ + O n can be obtained either by Schwarz 
reflection principle or by a similar calculation. We find 



L 



'E-+C, 
-Qn 



e 3 



^ ri _j^ n {^x, z)c dz 
" C Ai(£) Ai'(a^) - w 2 Ai'(f ) Ai(wry) 



e^dr/ + 0(n -1/<J e 



l/3 p -§|?| 3 ^ 



(234) 
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For the integral over (e + C n , oo) in (226), we need asymptotics of K n ^j >n (x, y). For x G E T ^/ 2 
and y € (e -+- C n , e + e), setting x = e + /3 _1 n -2 / 3 ^ and y = e + /3 _1 n~ 2 / 3 ry, we have 



(235) 



This follows from the analysis similar to that of (231). A weaker estimate is in (341), which is 



actually enough for our purpose. Hence for x 6 -Et,<e/2> 

K n . jjn (x,y)e n ^- 1 2 v ^dy 



e+C„ 



-Qn 



C 



£ - 7/ 



l/3^||5| 3/2 ^ 



(236) 



Now consider the integral over (e + e, oo). By the estimate (339) of K n -j jn {x, y) for x S Et )6 /2 an d 
y G (e + e, oo) and the identity 



we obtain 



V/n, !, n G(y;a) + H(y; Q )-2G(e;y) 

-V(y)+ay- -logQ n = . 

In I 



K n ^ n (x,y)e n ^- 1 2 v ^dy 
Qn jH O (n^e-f l«l 3/2 (l + M) -V(G( y ;a)-G(e;a))) dy = Q n0(e -f ICI 



(237) 



where 



e = min G(e; a) — G(y; a) > 0. 

j/>e+e 



(238) 



(239) 



Combining (236) and (238), we have 



e+C„ 



K n ^ n (x,y)e n ^ v ^di 



Qr 



C 



Ai(QAi / (r ? )-Ai / (g)Ai(r;) 



(240) 



We now insert the results (233), (234) and (240) into (226) and evaluate \fj n -j(x) for x E Er,e- 
The combination of the three integrals can be simplified. For this purpose, note that 



Ai^Ai'^-Ai'^AiQ?) 
£-77 

co 2 Ai(Q AiVr?) - Ai'(fl Ai^r/) 
£-77 

u Ai(Q Ai'^r?) - Ai'(g) Ai(o;77) 
£-77 



Ai(£ + t) Ai(r] + t)dt, 



Ai(£ + i) Ai(w 2 (77 + t))dt, 



Ai(£ + i)Ai(w(r? + t))eft, 



(241) 
(242) 
(243) 
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where we require "Str/ < !R£ in (242) and (242). This can be verified by using Ai (z) = z Ai(z) 
and the asymptotics [TJ 10.4.59 and 10.4.61] of Ai(z) and Ai (z) as z — > oo. Using the above Airy 
function identities, we find that 



AiCOAi'faJ-Ai'COAifa) 



c 



e ar >dr) 



oo roo 



Ai(£ + t) Ai(rj + t)dte ari dr] 



c Jo 

oo 



Ai(£ + t)e 



-at 



Ai(r/)e a ^ dt. 



Similarly, 

u? Ai(£) Ai'^) - Ai'(C) Ai(u 2 r]) 



c 



and 



C 



Ai(g) Ai' (uri) - u 2 Ai'(g) Ai(ur]) 



e ari diq 



Ai(£ + t)e 



c+t 



-at 



(244) 



c+t 



e ar >dri = u 2 / Ai(£ + t)e 



-of 



c+t 



Ai(u? f))e an dfi ) (it 

(245) 

Ai(ujf])e a ^dfi ] (it. 

(246) 



From these results and (|226j), (|233j), (|234j) and (|240j), we find that for x G £ T ,e/2, 



r n _,(a)^ n _,(x) = e «(— V(-)/2) _ Qr 



- w 



Ai(f + t)e~ at 
A\(u] 2 f\)e aTl dr) + w 



Ai(fj)e af ^dfi 



c+t 



C+t 
; 2 -o< 

C+t 
a 3 /3 



Ai{ujf))e an df) ]dt + o(l) 



. (247) 



Now the sum of three integrals inside the parentheses equals e a ' (cf. (223)), for all t. In order to 
see this, first note that the sum is independent of t since its derivative with respect to t equals 
from the Airy function identity 

Ai(z) +ojA\(ujz) +uj 2 Ai(u?z) = 0, z e C. (248) 

Then set t = and call the sum S(a). Taking the derivative of S(a) with respect to a and 
using (248) and the differential equation for the Airy function, we find that S'(a) = a 2 S(a). Now 
by noting that S(0) = 1 since Ai(fj)dfj = 1/3, we obtain that S(a) = e a3//3 . Hence 



T ri ^j(a)i) n -j{x) = Q r 



1 



e ax-\V(x) _ e a 3 /3 



uniformly for x £ £r,e/2- Therefore using (197) we find that 



1pn-j(x) 



#J>( e ) 



1 

Qn 



c -a 3 /3 1 c ax-hV(x) 



C_ a (0 + e- 



-a 3 /3 



Ai(C + t)e" Qt dt + o(l) 

Ai(£ + t)e~ a *(ft + o(l) 
<- s I +o(l) 



(249) 



— e ax-\V(x) _ p a£ 



(250) 
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uniformly for x = e + j3 l n 2 / 3 £ G -Et,<e/2- I n t ne l as t line, we used the identity 

roc 

/ Ai(£ + t)e- at dt = C- a (0- (251) 



e -« 3 /3 e a« 



□ 



Proof of (225). Let x > e + e/2. Using (338) and (214), a straightforward estimate implies that 



-y /-e+e/4 



n Je+Cn 



K n ^ n {x,y)e n{ay -^ V{y)) dy = 0(n-^ 2 (l + \ x \)~i e n(G(x;a)-H(x;a))/2y 



Similarly, we obtain using (336) 

roc 



1 

Qn 



e+e/4 



K n _ jin (x,y)e n(ay -^ V(y)) dy = 0((1 + \ x \yj e n(G(x;a)~R(x;a))/2y 



(252) 



(253) 



For the integral on T,± + C n , the calculation is easier than the proof of (224) since \x — z\ > e/2. 

(254) 



Straightforward estimates using Proposition 6.1 imply that 
1 



Qn 



T,±+C n 



K n _ iin (x, z)e naz dz = Oin- 1 / 2 ^ + \x\) 



-j n(G(x;a)-H(x;a))/2 



)■ 



Therefore, from (226), 



= e n(G(x;a)+K(x;a)-2G(e;a))/2 A + + ^^-^-ntH^-G^a))^ (255) 

= e n(G(x;a)+H( !r ;a)-2G(e;a))/2/ 1 + 



Hence using (197), we obtain (225). 



□ 



5.1.3 Proof 



Recall the outline (b2) in Section [2| 



From Lemma 
obtain 



5.2 



and Proposition 6.1l and using (12141) to estimate 7 ^ r e n ^ ax v ^/ 2 ^ — e a ^, we 



Wn-j,ll>n-j)lT 



C_ a (0 Ai(Ode. 



(256) 



T 



Now we evaluate := -7^{Kn-j,nXiT4>n-j)(e + P 1 n 2//3 £) defined in (78) asymptotically 



in L ([T, oo)). Using Lemma 5.2 and then estimating as in Subsubsection 5.1.2, we obtain 



K n -^ n {x,y)^ n -jdy = I!^f \ 



(K AhyX[ T,oo)C- a )(0 + 0(n-V3 e -§KI 3/2 ) (257) 
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for x = e + /3n 2 / 3 £ £ #T,e/2- For z > e + e/2, 



-j,n{x,y)ip n -jdy — 0{{1 + \x\) 



-j n(G(x;a)-H(x;a))/2 



The calculation is similar to Subsubsection 5.1.2 and we skip the details. Thus 



— r^r(-^Air y X[r,oo)C-a) 

/,nl e J 



0, in L 2 [T, oo). 



(258) 



(259) 



Therefore, from (79), 

P Jt _ J+1 , n (a;/J) = F (T) • (1 - (C_ a ,Ai) [T>oo) 



- ((1 - X[T,oo)-^AiryX[T,oo)) 1 -^AiryX[T,oo)C-a , Ai ) [T,oo)) 

Fi(T;-a). 



(260) 



Hence Theorem l.l [b) and Theorem 1.3^a) are proved. 
5.2 Proof of Theorem 11.51 

Note that G(xo(a c );a c ) = G(e;a c ) > G(x;a c ) for all x G (e, oo) \ {xo(a c )}. For a given in 



either (54) or (|56[), let xo( a ) be the point near xo(a c ) such that G(x; a) achieves its local maximum. 



The point xq(o) is well defined as long as \a — a c | is small enough. Note that for a > a c , xo(a) 



is same as in the definition of xq(cl) in Lemma 1.3 However, for a < a c , xq(o) is not defined in 



Lemma 1.3 We extend the definition of xq{cl) here for a < a c when a — a c is small enough. 



5.2.1 Proof of Theorem |1.5j|(a) 



We consider the double scaling situation 



a = a c + 



/3a 
^73 



■n 



(261) 



where a is in a compact subset of (— oo,0). Since we assume G"(xo(a c ); a c ) / 0, we have xq(o) 
^o( a c) = 0(\a — a c |) = 0(n -1 / 3 ). We also have, as in (97), using J^G(x;a) = x, 



G(x (a); a) = G(x (a c ); a c ) + x (a c )(a - a c ) + 0(\a - a c | 
= G(x (a c ); a c ) + ^^x (a c ) + 0(\a - a c | 2 ). 



(262) 



Hence since G(e; a) = G(e; a c ) + (a — a c )e by the definition of G, 

(3a 



G (xo (at); a) - G(e; a) = — -^(xq^slc) - e) + 0(n 



-2/3> 



We first evaluate r n _j(a) as in Lemma 5.1 Note that in Subsection 



5.1.1 



(263) 



we used properties 



of H(z;a) for the integrals over E± and properties of G(x;a) for the integral over (e, 00). Since 
there is no change in the properties of H(z; a), the integrals over S + and £_ are computed exactly 



47 



the same as given by (219) and (220). For the integral over (e, oo), note that the main contribution 
to (221 ) was from the part of y near e since G(y; a) takes its maximum for y near e. However, now 
due to (263) we need to add a contribution from y near xo(a). By using the standard Laplace's 

(264) 



method as in ( 99 ) , the contribution to the integral near xq (a) equals 

-^ > (x (a c ))e n ( G ^( a ) ;a )-^(l + o(l)) 



2tt 



Adding (264) and (221), 



nG"(x (a c );a c ))' 

X iin (x (a c ))e n(G(:Co(a);a) - G(e;a)) + o(l) ) . (265) 



2tt 



-G"(x (a c );a c 



But (263) implies that 

n(G(xo(a);a) - G(e; a)) = n 2 / 3 Pa(x (i 
for all large enough n since a < 0. Hence we find that 

Q 



e) + 0(n 1 / 3 ) < 



poo 

B j>n (e) / Ai(0e°*d£ + o(i; 
J o 



as in Lemma 5.1 Adding the integrals on the contours £±, we obtain 

Qr, 



T n -i(a) 



e^/ 8 (fl i)B (e) + o(l)), 



which is same as (197). 



The evaluation of ip n -j{x) is similar. We use the formula (226 ) as in Lemma 



5.2 



(266) 
(267) 

(268) 

The evaluation 



of the integral over (S + + C n ) U (£_ + C n ) is exactly the same as in Subsubsection 5.1.2 For the 
evaluation of the integral over (e + C ra ,oo) when x G E T ^/ 2 , we find that (236) is unchanged. For 

(269) 



the integral over (e + e, oo), the contribution near y = xo(a) implies that (238) becomes 

= n(G(:ro(a);a)— G(e;a)) 



K n _ j>n ( X ,y)e n ( a v-lnv)) d y = Q n O(„i/6 e -f l?l 3/2 1 



e+e 



n 



However, due to (266), this is again as in (p38|) . Therefore the result (12241) 

still holds for x G E Te / 2 . For a; £ (e + e/2,oo), the estimates (254) and (252) hold without any 
change. Moreover, it is straightforward to check that (253) still holds. Therefore (225) holds for 
x > e + e/2. Therefore, Lemma 5.2 holds without any changes. 



Now we proceed as in Subsection 



5.1.3 



In evaluating {ifi n -j, ^n-j)iT, the integral over (e + e/2) 



becomes, due to the contribution near xq(o), 



ji Yn—j I (e+e,oo) 



poo 

j 0(V^e n( - G{x ' a) - G( - e ' a)) )dx = o(e n( - G( - Xt,< - a ^ a) - G ^ a)) ) -> 

J e+e 



(270) 



by (266). This implies that (256) holds without a change. Similarly, it is straightforward to check 
that (259) holds. Therefore, we obtain (260) and Theorem l.E 'a) is proved. 
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5.2.2 Proof of Theorem |1.5|Kb) 

We consider the double scaling situation 



a = a c + 



a 



n 



where a' is in a compact subset of M. In this case, (262) and (263) are changed to 



and 



ct 

G(x (a);a) = G(x (a c );a c ) + x (a c ) h 0(n 

n 



Ot 

G(xo(a); a) — G(e; a) = — (xo(a c ) — e) + 0(n 

n 



(271) 



(272) 



(273) 



First we consider T n -j(a). There are two changes from the previous subsubsection. The first is 
that since a defined in ( |261 ) and a' defined in (271 ) are related as a = f3~ 1 n~ 2 ^ 3 a' , we have a — > 



and hence in (219), (220) and (265), we have a = in the integrals involving the Airy function. 



The second is that (|267|) does not follows from (265) since (266) no longer holds. Instead, due 

Qn 



to (|273|), (|265|) implies that 

2vr 



-G"(x (a c );a c )) 



. ^n(e) / AifcK 
Pv n \ Jo 

M jtn (x (a c ))e a '^ a ^ + o(l) 



(274) 



Hence adding (219) and (220) (with a = 0), we obtain 



r n -j(a) 



+ o(l) 



where 



2n 



-G"(x (a c );a c 



-A^ J> (x (a c ))e' 



a'(x (a c )-e) 



(275) 



(276) 



Now consider ip n _j(x). Like in Subsubsection 



5.2.1 



most of the estimates of Subsubsection 



5.1.2 remain the same. The only changes are the contribution near xo(& c ). Substituting (273) into 



(|269|) with a = (3 1 n 2 / 3 a', we have 

K n ^ n (x,y)e^y-^dy = Q n O( n -^e-^ m ) 



(277) 



for x £ E Te / 2 . This is similar to (238) except that e 2 n is replaced by n x / 3 . However, it is easy 



to check that this bound is enough for the rest of the analysis. The rest of the analysis of Subsub- 
section 



5.1.2 



continues without changes and we obtain the asymptotics (249) of T n -j(a)ip n -j(x 



for x G Et i€ /2 an d (255) for x > e + e/2. Using (275), we find that 



1pn-j{x) = Pj,n(a') 1 V '' 



Co(0 + e" a3/3 f-Le ""^ (*) _ e «A + (i) 



(278) 



49 



for x G Et i€ /2 an d 



^•(x) = Pj , n ( a ')^L e «(G(x;a) + H(,;a)-2G(e;a))/2 (1 + 
#j',n( e ) 



(279) 



for x > e + e/2, where 



^0 



D + D x (a') 



(280) 



The formulas (278) and (279) are different from (224) and (225) only by the factor pj jn (a'). 
We now prove the theorem. First, consider (58). From (279) and (330), we obtain 



(l/)n-j, lpn-j) jT( X0 ( ac )) =Pj,n( a ') 



■ f M iin (x)e n(G(x;a) - G(e;a)) (l + o(l))(i3 



BjA e ) J.jT(x {a c j) 

e n(G(a;o(o);o)-G(e;o)) roc 

=Pj,n{o ! , )-- 1 -rM j}n (xo(a c )) ^ =w= / 

&j,n{ e ) \f-nG "(xo(a);a)/27r Jt 

=(l-p j>n (a'))(l-G(T))(l + o(l)), 



e-3« 2 ^(l + o(l)) 
(281) 



by using ([272]) and g73| >. Also, for x € jj(x (a c )), by ( [336] ) and ( [27§ 



-j>Xjj(x (a c ))' i / ; n-i)( :r ) — / K n _j^ n (x,y)-ifj n -j(y)dy 

=0 ( e n(G(x;a)-M(x;a))/2j f Q^ e n(G(y;a)-G( e ;a))^ dy (282) 

= Q( e n(G(x;a)-H{x;a))/2^ 

Hence K n _j jn XjJ(z (a c ))' i /'™-.j ~~ ^ in L 2 ( jJ(xo(a c ))). Therefore, we find from (76) that 

Pn-i+l,n(«; (^o(a c ))) = 1 - (1 - p J >(o'))(l - G(T))(1 + o(l)) (283) 



and (58) is proved. 



Second, consider (57). We proceed as in the outline ["(b2)| in Section [2] We first evaluate 



(ipn-j,^n-j)iT = (ip n -j,if) n -j) ETie/2 + (ip n -j,ip n -j) ( e +e/2,oo)- For second term, a computation 



as 



in (281) yields that 



<Vn-;,VW-i>(e+ e /2,oo) = (1 ~ PjA^M 1 + ( 284 ) 

The first term is calculated as in Subsubsection |5.2.1| with the only change that the prefactor 
Pj,n{oi) is multiplied: 



j/ET,e/2 ~ Pj 



POO 

n (a0 J t c (6Ai(eK(i + o(i))- 



(285) 



Hence 



$n-j,i>n- j ) I T=Pj,n{a') / C (f ) Ai(£)<#; + (1 ~ Pj,n( a ')) + (286) 

JT 
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Now we evaluate := n X ' 2 Kn-j,nXiT'<pn-j{%) , x = e + /3 1 n 2//3 £ G /J. But it is straight- 

forward to check that, as in Subsubsection 5.2, 1[ 

Uj, n -PjA a ')l3 — r^(^AiryX[T,oo)Co) ^ in L 2 ([T, oo)). (287) 

Therefore, from ( |79| ), 

P n _ i+ljn (a;/J) =F (T) • (l - ( Pj>n (a')(C , Ai) [Tj0o) + 1 - Pi ,„(a')) 

-Pi,n(a')(( 1 - X[T,oo)^AiryX[r,oo)) _1 (^AiryX[T,oo)Co), Ai )[T,oo)) + o(l) (288) 

=^ > (a / )^i(r;0) + o(l). 



Thus (57) is proved. 



6 Summary of asymptotics of orthogonal polynomials and the 
Christoffel-Darboux kernel 

Define the matrix-valued function 

Y k (z;n)=( ^(^Pki^n) % 1 (p)(C<p k )(z; n) \ 

y-2vri7 fc _i(n)p fc _i(z;n) -2TTi-f k _ 1 (n)(Cif k ^i)(z;n) J 1 

for z G C\R. Here (Cf k )(z;n) is the Cauchy transform of (p k (x;n) = p k (x; n)e~ nV<yX \ The matrix 
Y k (z;n) is the solution to the Riemann-Hilbert problem for orthogonal polynomials with weight 
e -nV(z) j- gee ^20j). We are interested in the asymptotics of Y n -j(z; n) when j = 0(1) and n — > oo. 
We indicate the changes from the analysis of [17J and state the results. A similar derivation for the 
discrete weight can be found in [3]. 



Fix S > small enough. Let (see Figure 10 ) 



A s :={zeC\\z-e\<6}, (290) 
B s :={z G C | > e and \z - e| > 5} U {z G C | < e and |9fc| > 5}, (291) 

where e is the rightmost end-point of the support of the equilibrium measure. Comparing with 
notations in [17], As is the circle D €>aN+1 with 5 corresponding to the radius e. in [T71 Figure 1.4]. 
As in [T71 Figure 1.4], As is divided into four regions I, II, III and IV. Let Sjj be the contour in 
[T71 Figure 4.9]. We assume that the boundary of As is a part of and Bs is outside of the 
lens-shaped regions, cf. [T71 Formula (4.116)]. 

Several notations from |17j are used in this section, and we summarize them in Table [TJ Other 
notations may be slightly different but should be clear. 

By following the procedure of [17j . we find asymptotics of Y. Noting the symmetry 



Y(z) = I Y(z) I " , (292) 
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Figure 10: As and B$. 



Notation 


Definition in [17J 


n 


denned in Formula (1.21) 


9 


denned in Formula (1.24) 


7 


denned in Formula (1.26) 


u 


denned in Formula (1.29) 


d 


denned in Formula (1.30) 


u + (oo) 


explained below Theorem 1.1 




denned in Formula (1.34) 


P 


denned in Formulas (1.38) -(1.40) 




mentioned between Formulas (1.40) and (1.41) 


£« 


shown in Figure 1.5 


.y(°°) 


defined in Formulas (1.104)-(1.107) 


G 


defined in Formula (3.44) 


Si? 


shown in Figure 4.9 




defined in Formula (4.108) 



Table 1: Notations taken from [T7] 



it is enough to consider z G C+. The asymptotic formulas of Y(z) are different in B$ and Ag. For 
z G B s n C+, we have (cf. pTT, Formula 4.116]) 

Y n ^(z;n) = e% ff *R j , n (z)M$\z)e n <*W-& ff *. (293) 

The remainder Rj^ n (z) solves the Riemann-Hilbert problem similar to the Riemann-Hilbert problem 
for R in |17} Subsection 4.6] (cf. |17l Formulas (4.106)— (4.108)]) with the jump matrix which has 
the same estimate as in vr shown in |17} Figure 1.4]. This implies that uniformly for all z G C\£# 
and n (cf. \17\ Formula 4.115]) 

R j , n {z) = I + 0(n- 1 ). (294) 
The outer parametrix M^\z) solves the Riemann-Hilbert problem (cf. |17| Formulas (4.24)- 
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(4.26)]) 



M^'(z) is analytic in C \ 



MjZ\z)=(l + o(± 



z~> 
z j 



oo, 



z G 



(295) 
(296) 

(297) 



Note that the dependence on j in the asymptotics as z — > oo. The solution of this Riemann-Hilbert 
problem can be solved as in |17} Lemma 4.3]). Setting 

0(w + (oo) + d) 



6(oo) := 

e(z) ■, 



0(w+(oo) -d)' 
7(z)- 7 (z)-l0Mz) + d) 



it is straightforward to check that for z G Bg n C+, 



Mjj ) (z)=diag 



2^=0 



(298) 
(299) 



x diag 



Q(oo)- j 0(u + (oo) + d) 



Q{oo) j 9(u + {oo) + d) 



9(u + (oo) -£Sl- (2j - l)d) 9(u + (oo) + + (2j + l)d) J 

0(u(z)-^i1-(2j-l)d) » , >_j 7 _ 7 -i fl(u(^)+^n+(2j-l)d) \ 
U l Z i 2 e<u(z)+d) ™\ Z ) -2i " 6(u(z)-d) 



6(u(z)+d) 
■ 7 _ 7 -i 8(u(z)-^n-(2j+l)d) 



2i 



e{u(z)-d) 



e(u(z)-d) 

f^r^-i -y+y- 1 e(u(z)+^n+(2j+i)d) 
KJ \ Z ) 2 e(u(z)+d) 



and for z G Bg H C_ , 



1^1=0 a l+l 




x diag 



e(oo)- J '0(u+(oo) + d) 



e(oo) J '0(u+(oo) + d) 



0(u + (oo) 

-j: 



2tt 



- (2j - l)d) 0(u+(oo) + ^ft + (2j + l)d) 



g("(^)+^"+(2j-l)d) e{u(z)-^n-{2j-l)d) 
U l Z i -2i f?(u(z)-d) ^V*J -2 



-2i e(u(z)-d) ^ > -2 e(u(z)+d) 

,- 7+7 -i 9(u(z)+^n+(2j+i)d) Q^y y-gzl eM^)-^n-(2j+ii./i 



0(u(z)+d) 



-2( 



e(M(z)-d) 



(300) 



(301) 



Note that both formulas (300) and (301 ) can be extended to z G Bg n K. 

The asymptotics (293) especially implies the asymptotics of 7n-j- Since (cf. [Ml Formulas (3.10) 
and (3.11)]) 

-1/2 



7n-j 



-2ni lim z- n+ J +l (i?i,n(^M 1 ( ~ ) (z))i 2 e- n ( g ( z ) +£ 

|«|— s>oo ' J ' 



(302) 



we have from (294) and (300) that 



7„_ 3 - = 7n-i(l + 0(n *)), 



(303) 
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where 



In- j 



* (E^o a i+i ~ b i) 9(u + (oo) + d)e(u+(oo) + ^fi + (2j - l)d) 
2 9(u + (oo) (2j - l)d)9(u+(oo) - d) 



-i 



L«=o a i 



1=0 



h 



2j 



e(oo) 



(304) 



See [13 Formulas (1.62) and (1.63)] for the cases j = 1 and j = 0. 

Now consider z € A$. Then the analysis of the local parametrix as in |17[ Section 4.3] implies 
that (cf. [13 (4.119)-(4.121)] ) 

Y n „j(z;n) = e^R hn {z){M^ n ) p (z)e<^-^ (305) 

for z is in regions I and IV in |17} Figure 1.4], 



Y n ^(z; n) = e 2 a3 i? jin (z)(M jin ) p (z) I ^ nG{z) 1 



for z is in regions II in \17\ Figure 1.4], and 

nl 

Y n .j(z;n) = e^ a *R jin {z){Mj, n ) p (z) 



1 0' 

_ e nG(z) y 



e i(g(*)-a) ff 3 



n(g(z)- 5)0-3 



(306) 



(307) 



for z is in regions III in |17l Figure 1.4]. Here the local parametrix (Mj jn ) p is given by (cf. |17| 
Formulas (4.75) and (4.76)]) 



(DO)/ , 1 



(M jin ) p (z) :=M^>(z) 
where $(z) denotes <& aN+1 



2i\l 1 / 
z) in [T7]. We note that by definition 
3n 



(*(z)) ff3 / 4 ?($(z)). 



<&(*) 



-(2g(z)-F(z)-£) 



2/3 



so that 



®(z) = (3n 2/3 (z - e)(l + 0(|z - e|)) as z -> e 



(308) 



(309) 



(310) 



with /3 defined in (|22j) (see [13 Equations (1.34), (1.35), (4.74)]). 

We now summarize the asymptotics the orthonormal polynomials and their Cauchy transfor- 
mations. For notational convenience, we denote for z £ B$ 



M j>n (z) :=%. J ef(Mff) u (z), 



3>* 



Mj, n (z) :=7n-,e^(M)7) 12 (z), 



and for z £ As 



Bj,n( Z ) 



vn-^^ef ((Mff)n(z) - i(M^) 12 (z)) *{z?l\ 
™^%_ 3 ef (-(^)nW - i(M^) 12 (z)) *(*)-*/* 



(311) 
(312) 

(313) 
(314) 
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Remark 6.1. The formulas (300) and (301) contain the iV-variable Theta function 9. If N = 0, 
i.e., the equilibrium measure is supported on one interval, then 9 = 1, and the expressions of Jn-j, 
Aij in (z), M.j tn (z), Bj tn (z) and Vj :n {z) are much simplified. In particular M.j, n (z), Aij, n (z), B^ n {z) 
and Dj tTl (z) do not depend on n when N = 0. For example, for z G C+ 



M jn {z) =J — 2 , 7 + 7 1 i 1 7 ! V, (315) 



(316) 

and 

B iin (e) = ^(ai - 6o)~ 1/4 /? 1/4 . (317) 

Proposition 6.1. There exists 5o > sac/i f/iai /or eac/t /feed 5 G (0,5oL £ae following holds as 
n — > 00 and jf = 0(1). 

(a,) For z G £5, 

^•(z; n) = ^We^W-^H/ 2 ), (318) 
where Mj n [z) is an analytic function in B$ and 

M j>n (z) = M j>n (z)(l + Ofn" 1 )) (319) 

uniformly in z and n. The function Mj >n satisfies that (i) Mj^ n (z) = 0(z~ :) ) uniformly in 
n as z — > 00, fwj in any compact subset K G B$, Mj >n (z), M'- n {z) and \/Mj, n {z) are 0(1) 
uniformly in n and z £ K, and (Hi) Aij n (x) > and Mj n {x) > for all real x > e. 

(b) For z G 5 5 , 

(CVn-i)(*; n) = M j>n (z)e n ^^ +i / 2 \ (320) 
where Mj n (z) is analytic in B$ \R, continuous up to the boundary, and 

M j>n (z) = M j>n (z)(l + 0(n- x )) (321) 

uniformly in z and n where .A/lj in (z) defined in ( |312j ) is analytic in B$. The function ATj n 
satisfies that (i) Mj n {z) = 0{z ] ) uniformly in n as z — )• 00, (ii) in any compact subset 
K G Bg, Mj tn (z), Mj n (z) and 1/Mj >n (z) are 0(1) uniformly in n and z G K \ R and (Hi) 
—iM.j,n{x) > for x > e. 

(c) For z G 

^j(z) = (n l '« Ai(Q(z))B jtn (z) + K< (<S>( z ))D hn {z)) e~^ z \ (322) 

where Bj^ n (z) and Dj^ n {z) are analytic functions in Ag and 

B j>n {z) = B hn {z)(l + 0(n- x )), D jtn (z) = V j>n {z){l + O^ 1 )) (323) 

uniformly in z and n. The functions Bj >n and Dj tU satisfy (i) Bj >n (z), Dj tn (z), Bj n (z), 
D'- n {z), 1/Bj tn (z) and \/Dj >n {z) are 0(1) uniformly in n and z G As and (ii) Bj tfl (x) > 
and Bj tn (x) > for x G (e — 5, e + 6). 
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(d) For z G As n C+, we /iai>e 

(CVn-iX*) = e W3 (n 1/6 Ai(w 2 *(*))B iin (z) + n~ l / 6 u; 2 Ai'(oj 2 $(z))D j<n (zj) e~^ z \ (324) 
and /or z £ A$ (~) <C_ , we /taue 
(C<p n -j)(z) = -e m / 3 (n^oj 2 Ai(u$(z))B jtn (z) + n" 1 / 6 Ai'(u<f> (z))D j>n (zj) e~^ z \ (325) 



where uj = e 27 ™/ 3 and Bj n (z) and Dj >n (z) are the same functions in (322). The formulas (324) 
and (325) hold up to the boundary z G A§ n M. 



Proof. By formulas (300) and (301) of (z), the properties of the theta function 9 and the 

definition of d, we have that for p,q = 1,2, the functions (AC-„ and (M^') pq (z)~ 1 are 

uniformly bounded for z in any compact subset K C B$ and the functions 



m^UW + (M%\ 2 (z)) {z-ef/\ 

.-(J^VC*) + (M^) p2 (z)) (z - e)~^\ 
(m^U(z) + (M^U^y 1 (z - e)^ 



are uniformly bounded for z G A$ and n. We note that although Mj°^\z) is not well defined 
on J, the functions considered above are well defined on Ag n J. Plugging in these estimates 
and the estimate (294) of Rj jTl (z) into (293), ( |305 ), (306) and ( 307| ) we obtain the estimates of 
Mj, n (z), 1/M j>n (z), Mj, n {z), 1/Mj, n (z), B jj7l (z), 1/B j>n (z), D jtn {z) and 1/D jj7l {z). With the help 
of the Cauchy's integral formula, we further obtain the estimates of M'^ n (z), M'^ n (z), Bj n (z) and 
Dj n (z). From the formula of M^\z) we also derive the positivity of Mj >n (x), Mj iH (x), —iMj in (x), 
Bj <n (x) and Bj >n (x). □ 

Remark 6.2. We use the following identity in the analysis. It is straightforward to derive from the 
Riemann-Hilbert problem of Y/~(z;n) that det Yfc(z;n) = 1. This implies that 



Ik 
Ik-i 



-1 

2Tri 



{Pk{z)C<Pk-i{z) -pk-i(z)C(p k (z)). 



(326) 



Taking k = n — j and using asymptotic formulas (322), (324) and (325) in (326), with the help of 
[H 10.4.11 and 10.4.12] 



we find 



M(z)u Ai' (uz) - Ai(z) Ai(uz) -- 
Ai(z)w 2 Ai'(w 2 z)-Ai'(z)Ai(w 2 z) = 

Bj >n (z)Dj + i 7 n(z) - Bj+i t n(z)Dj >n (z) 



2n ' 

= 7ri/6 

ln-j 
ln-j-1 



(327) 
(328) 

(329) 
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Proposition 6.1 implies the following asymptotic properties of ip n -j- These are used in the main 
analysis extensively. 

Corollary 6.1. Fix T £ R. There exists 5q > such that for each fixed e G (0,5o]; the following 
holds as n — > oo and j = O(l). 



(a) For x > e + e, 



^n-j(a:) =M J> (x)e n ( G W- H ( x »/ 2 = 0(e n < G ( s )- H ^)/ 2 (l + 



forM jtTl (x) in ( |318| >. 

f&j Ze£ i?T,e := [e + /3 -1 n~ 2 / 3 T,e + e] 6e #ie interval defined in (116). For x G -Er, £ , 

^n-i(x) = O^^e-f Kl 3/2 ), £ := /3n 2 / 3 (x - e). 

(c) Let 1^ := [e + f3~ 1 n~ 2 ^ 3 T, oo) 6e i/ie interval defined in (23). JTien 

llV'n- i ||L 2( /T ) = 0(n- 1 / 6 ). 
Afeo /or every x in (e + e, oo), i/iere is e' > suc/i that 

ll^n-illia([s,oo))= 0(e" e ' n ). 
('dj ^4s n — > oo, Vj,n(0 := "0n-i( e + /3 _1 n _2//3 ^) satisfies 



(330) 



(0-^Bi,n(e)Ai(0, m L 2 ([T, oo)). 



(331) 



(332) 



(333) 



(334) 



Proof. 



G(z; a) — H(z; a). 



he result (a) follows from (318) and noting that 2g(z) — V(z) — 

(b)l note that T < £ < en 2 / 3 . Thus, |Ai(f)| < Ce-§^ 3/2 and \Ai'(Q\ < C(|C| 1/4 + 



For 

l)e-§IC~< C'nVee-fl?! 372 for some constants C, C > 0. From (|322l) and the behavior of <3? 



in A e , we obtain the estimate with the factor | changed to a smaller constant which can be made 
arbitrarily close to I if we take e smaller. To be definite, we fix this constant as I. 



For (c) by the asymptotics (330) and (331) of tp n -i, we find 



<2 



-j\\ 2 L 2 (E T , t ) + ll^n-illi 2 ([e+e,oo)) 



/3n 2 / 3 < 



IT 

--Oin- 1 / 3 ). 



Pn 2 / A J e+t \ ) 



The estimate ([333) is similar. 

Item |(d)] follows from Proposition 6.1 (c). 



(335) 



□ 
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The above asymptotics for ip n -j yields the asymptotics for the Christoffel-Darboux kernel 
K n ..j >ri (x,y). 

Corollary 6.2. Let T G R be fixed. There exists 5q > such that for each fixed e G (0, <5q], the 



fallowings hold as n — > oo and j = 0(1). £e£ -Et,e &e £/ie interval defined in (116). 
(aj For x,y £ (e + e/2, oo), 



ir n _ iin (x,2/) = ( e n(G(*)-H(x)+G(y)-H( W ))/2 (1 + j^-;^ + 



y\ 



(b) Forx,ye E T>€! 

K„_ Jin (x, y ) = 0(n 2 / 3 e-f(l«l 3/2 +H 3/2 )), 
where £ := /3n 2//3 (x — e) and 77 := /3n 2 / 3 (y — e). 
(cj For 3; G (e + e, 00) and y G E Te i 2 , 

K n _ j>n (x,y) = 0( n l/6 e n(G(x)-H(,))/2 e -||r,|3/ 2(1 + j^-^ 

(c2) For x G E T ^i2 and y G (e + e, 00), 

y) = 0(n 1 /6 e -|l5l 3 / 2 e n(G(,)-H(,))/ 2(1 + |y|ri) _ 

All estimates above are uniform in x, y in their domains and in n. 



(336) 



(337) 



(338) 



(339) 



Proof. Item (c) and 1(d) I follow directly the asymptotics (330) and (331 ) of iftn-j, and the Christoffel- 



Darboux formula ( |69[ ) of K n _j jn (x,y), noting that x — y never vanishes. 
For x, y G (e + e/2, 00), ( p30| ) impl ies that 



TC U oA - MG(x)-H(x)+G(y)-H(y))/2 M jA X ) M j+lAy) - M j+l,n(x)Mj, n (y) 

n— j ,n\X , y) — c 



x-y 



(340) 



Since Mj :U and its derivatives are uniformly bounded, we obtain (a 



Item |(b)] follows from a similar calculation but using the asymptotics (322). The calculation is 
direct and is the same as |14} Formula (3.8)]. □ 

We also need the following results for the Christoffel-Darboux kernel. 

Corollary 6.3. Fix T G M. and let 1% be the interval defined in (23). Then we have the following: 

(a) For any fixed C, we have 



/3n 2 / 3 



K n _ j>n {x, y) = K AiTy (£, 77) + o(e 



(341) 



for all x,y G 1^ as n —> 00, where £ := (x — e)/3n 2 / 3 and n := (y — e)/3n 2 / 3 . 



58 



(b ) Define the operator /C n _,- n by kernel 

ICn^n&rj) := ^K n ^ n (e + + -^_). (342) 

Then 

j,nX[T,oo )) (1 -X[T,oo)^AiryX[T,oo)) (343) 

in trace norm as n — > oo. 

(qj TTie operator norms of (l — Xi T ^n-j,nXl T ) 1 are bounded uniformly in n. As a corollary, 
The operator norms of 

(1 ~ X[i n ,oo)-^n-i,nX[j; n ,oo)) , (344) 

are a/so bounded uniformly in n and in x n as long as x n are in a compact subset of (e, oo). 
fay We have 

lim det (l - XirK n -j, n XrA = *b(T), (345) 

and 

n^o det i 1 ~ X[x,oo)K n -j,nX[x,oo)) = 1 - (346) 
for any x is in a compact subset of (e, oo). 

Proof The proof of a result similar to | (a) | for the non- varying weight is given in |14t Formula (3.8)] 
. The varying weight case is proved in the same way. Note that our C is the c in [14, Formula 
(3.8)], which can be assumed to be an arbitrarily large positive number. 

The proof of a result similar to |(b)| for the non- varying weight is given in the proof of the (5 = 2 
case in |14| Corollary 1.4] . The varying weight case is proved in the same way. 

Items |(c)| and |(d)| follow from |(b)[ 

□ 
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